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I.  INTRODUCTION 

The  contract  entitled,  "Investigation  Of  Electromagnetic  Coupling 
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75-C-0293)  was  awarded  to  the  University  of  Illinois  for  one  year,  starting 
December  1,  1977.  The  contract  was  later  extended  to  Feb.  28,  1979  at  no 
additional  cost.  This  is  the  final  report  of  the  contract. 
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ABSTRACT 

This  paper  addresses  the  problem  of  computing  the  current  induced 
in  a thin  wire  located  inside  a cylindrical  cavity  with  a circumferential 
slot  when  the  cavity  is  illuminated  by  an  incident  plane  wave.  The  calcu- 
lation  is  carried  out  in  two  steps.  First,  the  problem  of  penetration 
of  the  incident  field  into  the  cavity  is  solved  by  the  method  of  moments 
under  the  assumption  that  the  presence  of  the  wire  inside  the  cavity  creates 
little  or  no  perturbation  of  the  interior  field.  Next,  the  induced  current 
on  the  wire  is  calculated  by  the  following  two  methods:  (i)  use  of  a simple 
analytical  formula  derived  from  the  application  of  the  Wiener-Hopf  techniques 
to  the  finite  wire  problem;  (ii)  numerical  solution  of  an  integral  equation. 
Extensive  numerical  results  for  the  induced  current  are  presented.  It  is 
found  that  the  current  is  sensitive  to  the  cylinder  radius,  the  cavity 
height,  the  frequency  of  excitation,  and  the  wire  location,  but  is  relatively 
less  sensitive  to  the  variation  in  the  slot  length.  In  addition,  the  induced 
current  on  a wire  inside  the  cavity  can  be  much  larger  than  its  counterpart 
in  free  space  illuminated  by  the  same  incident  plane  wave  at  frequencies 
where  the  cavity  is  near  resonance. 
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I . INTRODUCTION 


The  penetration  of  an  electromagnetic  wave  through  an  aperture  into 
a cylindrical  structure  is  of  current  interest  because  of  its  application 
to  EMP,  to  EMC,  and  to  biological  studies.  As  early  as  1949,  Sotmnerfeld 
[1]  studied  the  problem  of  an  infinitely  long  circular  cylinder  with  a 
longitudinal  slot,  illuminated  by  a normally  incident  plane  wave.  Using  a 
Fourier  analysis  approach,  Sommer feld  reduced  the  problem  to  a system  of 
infinitely  many  linear  equations,  but  declared  resignedly,  "We  can  do 
practically  nothing  with  the  problem."  Silver  and  Saunders  [2]  used  the 
saddle-point  integration  method  for  the  inversion  of  Fourier  transforms 
and  extracted  the  far  field  of  Sommerfeld's  problem.  Hitherto,  a number 
of  extensions  along  this  line  of  work  have  been  reported  [3],  [4].  With 
the  advent  of  high-speed  digl.-al  computers,  Sommerfeld's  penetration 
problem  can  now  be  solved  by  numerical  means  [5],  [6].  The  penetration  of 
an  EM  wave  into  the  cylinder  through  a rectangular  aperture  was  first 
carried  out  by  Safavi-Naini,  Lee  and  Mittra  [7],  [8].  Their  problem  has 
a more  complex  geometry  than  Sommerfeld's  in  that  two  conducting  plates 
are  introduced  inside  the  cylinder  at  z = =h  to  form  a cavity,  as  illus- 
trated in  Figure  1. 

In  the  present  report,  and  extension  of  the  penetration  problem  of 
Safavi-Naini  £t  al.  is  investigated;  in  which  a thin  wire  is  added  inside 
the  cavity  and  the  problem  is  to  determine  the  current  induced  on  the  wire. 
The  wire  is  oriented  parallel  to  the  longitudinal  direction  of  the 
cylindrical  cavity.  If  the  slot  in  the  cavity  wall  is  also  longitudinal, 
there  is  little  induced  current  on  the  wire.  Hence,  we  concentrate  on 
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the  more  interesting  case,  namely,  the  slot  which  is  circumferential  on  the 
cylindrical  cavity  wall.  The  composite  geometry  of  the  present  problem 
is  sketched  in  Figure  2. 

Due  to  the  thinness  of  the  wire,  it  appears  reasonable  to  assume  that 
the  presence  of  the  wire  does  not  perturb  the  field  generated  Inside  the 
cavity.  Thus,  the  problem  under  consideration  can  be  solved  in  two  steps: 

(A)  Determine  the  field  E inside  the  cavity  as  if  the  wire  were 
absent. 

(B)  Using  E as  an  incident  field,  determine  the  Induced  current 
I(z)  as  if  the  wire  were  situated  in  the  free  space. 

We  emphasize  that  the  above  two-step  approach  is  an  approximation.  The 
exact  degree  of  approximation  will  be  studied  in  a separate  report. 

The  plan  for  the  present  report  is  as  follows:  In  Section  II,  the 
problem  of  Part  (A)  with  the  wire  absent  is  formulated  and  a system  of 
infinitely  many  linear  equations  derived.  The  procedures  are  briefly 
described  below: 

1.  The  unknown  electric  field  across  the  aperture  is  represented 
by  a Fourier  series. 

2.  Applying  the  equivalence  principle,  the  aperture  is  shorted  by 
a perfect  conductor.  The  effect  of  the  original  aperture  field 
is  accounted  for  by  introducing  equivalent  magnetic  currents  on 
both  sides  of  the  shorted  aperture. 

3.  Inside  the  cavity,  the  magnetic  field  produced  by  the  equivalent 
magnetic  current  is  determined  via  a magnetic  vector  potential. 

The  resultant  field  is  given  in  the  form  of  a doubly  infinite 
series  of  eigenfunctions  of  the  cavity. 
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4.  Exterior  to  the  cavity,  the  magnetic  field  is  also  generated 
via  a magnetic  vector  potential,  which  is  partitioned  into  three 
components  to  represent  the  incident  field,  the  reflected  field, 
and  that  produced  by  the  equivalent  magnetic  current. 

5.  The  coupled  equation  from  which  the  unknown  aperture  electric 
field  can  be  determined  is  developed  by  enforcing  the  continuity 
of  the  tangential  magnetic  field  across  the  aperture. 

6.  The  coupled  equation  is  solved  numerically  by  the  method  of 
moments  [9]. 

In  Section  III,  the  current  induced  on  the  wire  inside  the  cavity 
when  it  is  illuminated  by  the  field  E computed  in  Part  (A)  is  derived. 

Since  E is  given  in  terms  of  a doubly  infinite  series  of  eigenfunctions 
of  the  cavity,  it  can  be  interpreted  as  a spectrum  of  plane  waves.  As  a 
result,  the  induced  current  I(z)  can  be  determined  by  superimposing  the 
currents  due  to  each  component  of  the  plane  wave  spectrum.  The  current 
induced  by  each  component  of  the  spectrum  is  determined  by  one  of  the 
following  two  methods: 

1.  The  standard  numerical  (moment)  method  based  on  an  integral 
equation  formulation  [10],  and 

2.  the  simple  approximation  formula  recently  developed  by  Chang, 

Lee,  and  Rispin  [11],  [12]. 

It  should  be  noted  that  both  of  the  above  methods  apply  only  if  the  incident 
plane  wave  is  homogeneous.  In  the  present  problem,  however,  the  field  E 
inside  the  cavity  consists  of  both  homogeneous  and  inhomogeneous  plane 
wave  components.  Hence,  the  above  two  methods  have  to  be  extended  by 
analytical  continuation  to  cover  the  case  of  an  incident  inhomogeneous 
plane  wave. 
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In  Section  IV,  possible  difficulties  of  numerical  computations  of 
the  equations  derived  in  Sections  II  and  III  are  considered.  Techniques 
for  improving  the  convergence  rate  of  the  suimnation  procedures  are 
presented  in  Section  IV.  Extensive  numerical  results  are  presented  in 

Section  V. 


II.  FIELD  IN  THE  CAVITY 


In  this  section,  the  field  E excited  in  the  cavity  with  the  wire  absent 
is  derived.  The  geometry  of  the  problem  under  consideration  in  this  section 
is  illustrated  in  Figure  2,  with  the  wire  removed.  The  conducting  circular 
cylinder  is  infinitely  long  and  is  of  radius  a.  A cylindrical  cavity  is 
formed  inside  the  cylinder  by  two  conducting  plates  located  at  z * ±h. 

The  cavity  is  coupled  to  exterior  excitations  through  a circumferential 
slot  on  its  cylindrical  wall.  The  rectangular  slot  is  of  dimension 
2c  x 2d  and  centered  at  (x  - a,  y - 0 , z ■ 0) . The  width  of  the  slot  is 
assumed  to  be  small  in  terms  of  wavelength,  i.e., 

2kd  « 1 . (2-1) 


The  structure  is  illuminated  by  a normally  incident  plane  wave  of  unit 
magnitude  described  by 


— i A ikx 

E • ze 


1 ikx 
y_e 


(2-2) 


where  the  time  harmonic  factor  exp(+jwt)  has  been  suppressed  and 

H ■ /y/e  ■ 120  TT  is  the  intrinsic  impedance  of  the  free  space. 

The  symmetry  of  the  configuration,  together  with  the  plane  wave 

excitation  given  in  Eq.  2-2,  dictate  that  the  tangential  electric  field 

across  the  aperture  be  an  even  function  of  0.  Furthermore,  the  narrowness 

of  the  slot  enables  us  to  assume  that  the  aperture  field  is  approximately 

constant  in  z and  is  z-directed.  Thus,  the  tangential  electric  field  across 

the  aperture  E ■ E z can  be  represented  by  a Fourier-cosine  series: 

4 z 

00 

E - l E cos  r ♦;  |*|  < * |z|  <d  . (2-3) 
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In  accordance  with  the  boundary  condition  that  vanishes  at  <p  ■ -$q($q  m c/a), 


O'  } are  found  to  be 


! y " 0 ,1 ,2 , . . , 


(2-4) 


Due  to  the  assumed  direction  of  the  aperture  electric  vector  (or  the 


equivalent  magnetic  current) , it  can  be  shorn  that  a field  TM  with  respect 
to  z is  sufficient  torepresent  the  total  field  inside  the  cavity.  Such 
a field  can  be  generated  via  a z-directed  magnetic  vector  potential 

A 

A » A z.  The  relations  between  A and  the  field  components  are 

2 Z 
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(2-5) 


Inside  the  cavity,  A is  a linear  combination  of  all  possible  elementary 

z 

wave  functions  that  are  solutions  to  the  scalar  Helmholtz  equation  [13]. 
It  assumes  the  form 


A ( P , <(> , z ) - H A J (Y  P)cos  mtf)  cos  a z 
z r,  ran  m n n 

m,n*0 


(2-6) 


where 


.2  _ 2 
k 1 an 


n T2  . 2 . 2 
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ot  s ■ m * 0 jl  1 2 j • i 
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(2-7) 


(2-8) 


The  constants  (A  } are  unknowns,  and  {J  } are  the  mt'1  order  Bessel 
ran  on 

functions  of  the  first  kind.  Note  that  with  the  choices  of  the 

eigenvalues  {a^}  in  Eq.  2-8,  the  boundary  conditions  that  Ep  and  E^ 

be  zero  at  z » ±h  are  automatically  satisfied. 

The  constants  {A  } in  Eq.  2-6  are  unknown  quantities,  and  by 
nn 

enforcing  the  conditions  that  E be  zero  on  the  cylindrical  wall  and 
equal  to  Efi  on  the  aperture,  they  can  be  evaluated  in  terms  of  {E^}, 
the  expansion  coefficients  of  Efl  given  in  Eq.  2-3.  First,  we  evaluate 
E (referring  to  Eq.  2-5)  at  p ■ a to  obtain 


w 

X 2 

y y y A J (y  a)  cos  m<f>  cos  a z 

laipe  Ln  n mn  m n " 

m,n*u 


G(<|>, z)  l E cos  rj 
U-0  M M 


|4>|  < it  , |z|  <_  h 


(2-9) 


where  Q is  the  characteristic  function  of  the  aperture: 

i ; f<M  i $0  . M 1 d » 

0 ; otherwise. 

We  recognize  that  Eq.  2-9  is  a standard  Fourier-Bessel  series  with 


n(<J>,z) 


(2.10) 


unresolved  coefficients  {A  }.  Hence,  {A  } are  determined  by  standard 

mn  mn 


procedures  with  the  results 


i4um£  cos  m<J>rt  sin  a d <» 

0 n r _ 

A * 5 ) E 

e G to  y hJ  (y  a)  y*0 
m n n'n  m n 


(-Duru 

r2-m2 

u 


(2-11) 


m,n  * 0)1)2) i 


in  which  e is  the  Neumann  number,  defined  by 
m 


[2,10-0 

m V 1 , m 0 
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In  the  region  exterior  to  the  cylindrical  structure,  the  symmetry 
of  the  configuration  again  leads  to  a TM  field  with  respect  to  z.  Thus, 
the  exterior  field  is  also  determined  via  a z-directed  magnetic  vector 
potential  A*  ■ ¥z.  For  reasons  which  will  be  obvious  later,  ¥ is 
partitioned  into  three  wave  funcitons:  ¥ » ¥*  + ¥r  + ¥8.  The  first  of 
these  ¥i  represents  the  incident  plane  wave;  it  is  independent  of  a and 
is  given  by 

00  Q 


¥*  ■ ^ 7 J (kp)cos  m<)> 

fit  *•  r m 


(j)  •*-  e m 

m»0  m 


(2-12) 


The  field  components  of  ¥*  are  given  in  Eq.  2-2.  The  second  wave  function 
¥r  represents  the  reflected  wave  when  the  aperture  is  closed  by  a perfect 
conductor;  it  assumes  the  form: 


.“.ml  J (ka) 

l H--* 


“ m^O  Em  \ H(2) (ka) 


H^(kp)cos  nxt> 
01 


(2-13) 


m 


where  {H^}  are  the  mth  order  Hankel  functions  of  the  second  kind.  The 
m 

third  wave  function  ¥S  corresponds  to  the  field  produced  by  the  equivalent 

^ r s 

magnetic  current.  In  contrast  to  both  ¥ and  ¥ , ¥ depends  on  a and  is 
represented  by  a continuous  spectrum  of  cylindrical  waves : 


¥s  - ? cos  1$  If  (a)H( 

m“0  i ” ” 

where  {F  } are  unknown  functions,  and 


(2)(yp)ejotz  da 


2 2 

k > a 

2 2 

k < a 


V + 

To  determine  (F  } , we  first  evaluate  E at  p * a to  obtain 
01  z 


(2-14) 


(2-15) 
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(2-16) 


jo cos  ^ r Y2va)Hi2)(Ya)ejaz<ia 


- !i(<p,z)  l E cos  r 4 , |<fr|  <_  IT  , |z|  _<  ® 

p-0  p M 

Then,  we  multiply  both  sides  of  the  equation  above  by  cos  n<J>,  n - 0,1,2,..., 

and  after  that,  integrate  the  equation  over  the  entire  domain  of  interest. 

By  invoking  the  orthogonal  properties  of  (cos  n<J>},  (F  } are  determined 

in 

in  terms  of  {E^}: 


j2u)ye  cos  m<t>  sin  ad  00  (-l)^r 

F-(a)  ” 2 2 (2)  £ Fu  2 2 

e 1Toy  Hu;(ya)  y-0  P T - m 


m 


(2-17) 


m 


m 


m—  0,1,2,...  . 

With  {A  } and  {F  } defined  by  Eqs.  2-11  and  2-17,  tne  requirement 
mn  m 

that  the  tangential  electric  field  be  continuous  across  the  aperture  is 

automatically  satisfied.  However,  these  definitions,  of  themselves,  do  not 

ensure  the  continuity  of  the  tangential  magnetic  field  across  the  aperture. 

To  enforce  the  continuity  of  the  magnetic  field,  we  proceed  as  follows:  In 

the  region  exterior  to  the  cylinder,  the  three  partial  magnetic  fields 
i r 3 

H^,  H^,  and  (Hp  is  immaterial)  corresponding  respectively  to  the  three 
wave  functions  defined  in  Eqs.  2-12  through  2-14  are 


H 


H 


H 


i 

4> 


r 

♦ 


s 

♦ 


«>  m 

1 j'(kp)cos  m4>  , 

joiy  Ln  e m Y ’ 

J m-0  m 


m 


-2k  y J I 
m-0  £m 


Jm(ka) 

H(2)(ka) 

in 


H^2\kp)cos  m<t> 

m 


— I 


cos 


m-0 


f yf  (<W 

J m m 


(2)'(Y0)ejaZ  da 


(2-18) 
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Inside  the  cavity  H is  given  by 


H - — - y y y A J (y  p)cos  m<f>  cos  <S  z . 
d>  u L Ln  n ran  m n n 

y m,n*0 

Hence,  the  desired  continuity  of  the  tangential  magnetic  field  across 
the  aperture  now  reads 


(2-19) 


oo  m 


2k.  y 1- 

* Ml  I *>  P 


Jwv 


mpO  em 


J_(ka)  ^2V 
J’(ka)  --=« H<2)  (ka) 

m H(2)(ka)  m 

m 


cos 


--  y cos  m<j>  yF  (a)H^  ' (ya)e^aZ  da 
u mm 


mpO 


00 

- — yy  yA  J'(y  a)  cos  m4>  cos  a z 
u n mn  m n n 


(2-20) 


m,n*0 


|4»|  1 1*1  i d • 

The  above  equation  can  be  simplified  by  recognizing  that  the  term  in  the 
bracket  on  the  LHS  is  the  Wronskian  of  Bessel  functions: 


J_(ka) 

J*  (ka)  - H^Z),(ka) 

m „(2),,  m 


J2_ 


Hv  ' (ka) 
m 


irkaH^  (ka) 
m 


By  substituting  the  above  result  into  Eq.  2-20,  it  becomes 


J cos  m4>  [ yF  (a)H^' 
u - J m m 

®*0  > 


(2)'(ya)e:,az  da 


y T cos  m<t>  y A J’(y  a)  cos  a 2 
L L.  n mn  m n n 

m,n»0 


4 i ; 1*1  i*a-  1*1  - d • 

wa  m-0  e Hu;(ka) 
in  m 


(2-21) 
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I ' 


■ ■ 

■ : 


Applying  the  method  of  Galerkin  [9],  the  preceding  equation  can  be  solved 
numerically.  First,  we  integrate  the  equation  with  respect  to  2 over  the 
aperture;  the  result  is 


l cos  m$  f*  2 3l-"-otd  yF  (a)H^2\ya)  da 
Ln  J a mm 


« 2 sin  a d 

- I l cos  m<J>  - — — Y_A  J' (Y  a) 

*•  *•_  an  mn  m n 

m,n“0  n 


m m 


(2-22) 


Next  we  multiply  both  sides  of  the  equation  above  by  cos  v<j>,  V ■ 0,1,2,..., 
and  the  integration  of  the  resultant  equation  over  the  entire  domain  of 


interest  leads  to 


l g_  f Slna-  YF  (a)H^2)'(ya)  da 
mV  J a m m 


“>  sin  a d 

- 7 7 g — y A J'  (y  a) 

L Ln  smv  a 'n  mn  m n 
m,n*0  n 


CD  m Q 

_4d  c °mV 

TTwa  e „(2)  . 

m-0  m H (ka) 
m 


; V « 0,1,2,. 


(2-23) 


where  g stands  for 
mv 


sin(m  - v)<I>q  sin(m  + v) <J>q 
®mv  ” m-v  + m + v 


(2-24) 


Then,  we  replace  {A  } and  {F  } in  Eq.  2-23  by  their  definitions  in  Eqs. 
mn  m 

2-11  and  2-17.  After  some  algebraic  manipulations,  we  arrive  at 


I E A^+  A^  « D ; v * 0,1,2,... 
li  I yv  pv  vj 


(2-25) 


we  intended  to  derive.  In  general,  this  system  of  equations  cannot  be 
solved.  However,  if  the  series  in  Eq.  2-3  is  truncated  at  a finite  number 
N,  the  system  would  become  a system  of  N x N algebraic  equations,  with 


14 


III.  CURRENT  INDUCED  ON  A WIRE  IN  THE  CAVITY 


In  this  section,  the  current  induced  on  the  wire  inside  the  cavity 
when  it  is  illuminated  by  the  field  E is  derived.  We  attack  the  problem 
by  assuming  that  the  presence  of  the  wire  does  not  perturb  the  field 
generated  inside  the  cavity.  It  enables  us  to  use  E as  an  incident  field 
and  to  determine  the  induced  current  I(z)  as  if  the  wire  were  situated 
in  the  free  space. 

With  reference  to  Eq.  2-26,  we  note  that  the  electric  field  tangen- 
tial to  the  wire  Ez  is  given  in  terms  of  a doubly  infinite  series  of 
eigenfunctions  of  the  cavity.  The  series  can  be  interpreted  as  a spectrum 
of  plane  waves.  Explicitly,  Ez  is  rewritten  as 

00 

E - 7 f (p,4>)  exp(+jk  cos  9 z)  ; (3-1) 

2 n n 

n»- oo 

p <_  a,  |(J>|  _<  ir,  | z | £ h , 

where 

en=C0S_1(kE|  * (3_2) 


and 


— I 

TTh  L 


(-1) Pr  sin  a d °°  cos  m<t>  cos 

E M 3-1 

U*0  ^ an  m-0 


t>0  J (YnP) 
U m n 


E (r2  - m2)  J (y  a) 
my  m n 


‘ (3-3) 


We  interpret  each  component  of  the  field  in  Eq.  3-1  as  a plane  wave  in 
free  space  propagating  in  the  direction  0^  with  respect  to  the  z-axis. 

In  this  report,  two  methods  are  used  to  compute  the  current  induced 
on  the  wire  due  to  each  component  of  the  plane  wave  spectrum  described 
in  Eq.  3-1.  The  first  of  these  methods  is  based  on  an  integral  equation 
formulation: 
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z^+h 


In  bn  ' II 

SH  i 

i Z_-h 
0 w 


)K(z,z')dz'  - - j4iro)cf  e-l^cos  ®nz  . (3-4) 


The  above  thin  wire  scattering  problem  has  been  thoroughly  studied  in 
recent  years,  and  a number  of  efficient  programs  to  compute  the  unknown 


• (1) 


current  1^  ' have  been  developed.  The  program  developed  by  Butler  [10], 


which  is  based  on  solving  Eq.  3-4  by  the  method  of  moments,  is  adopted 
here.  Applying  the  principle  of  superposition,  the  induced  current  I 
is  given  by 


I(z)  - l ffp.OOi^U)  , 

u n n 

n“-°° 


(3-5) 


where  p is  evaluated  at  the  location  of  the  wire. 

As  is  well-known,  using  moment  methods  to  compute  the  current  on  a 
wire  is  extremely  time-consuming  when  the  wire  is  of  the  order  of  several 
wavelengths.  An  alternative  method  that  is  suitable  for  long  wires  is 
to  use  the  simple  approximation  formula  developed  recently  by  Chang,  Lee 
and  Rispin  [11],  [12].  The  techniques  used  to  derive  this  simple  formula 
are  briefly  described  below: 

1.  Using  a Wiener-Hopf  method,  the  reflection  coefficient  from  the 
end  of  a semi-infinite  wire  illuminated  by  a plane  wave  of  unit 
amplitude  is  determined. 

2.  By  considering  the  multiple  bounces  of  the  current  waves,  the 
current  Induced  on  a wire  of  finite  length  can  be  expressed  in 
terms  of  two  Neumann  series  involving  the  just  mentioned  reflection 
coefficient.  The  series  are  then  summed  up  into  a closed  form 

to  give  the  desired  approximation  formula. 
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*•  * 


The  induced  current  due  Co  each  component  of  Che  plane-wave  spectrum  is 

.(2) 


denoted  by  I and  is  given  by 
n 


(2)  f"  ~z') 

[,  ROT.Zh,)  ^(^9n»hw-Z')V(9n»hwjU  V> 

r ROr,h+*')  “] 

+ C(6  ) — „/  „ \ + A(0  ,h  +*•>  V (tt-0  ,h  ) U (h  +z ' ) 

I n R(ir,2hw)  n w n w\  w 


+ V(0n,z*) 


(3-6) 


where  z'  • z - Zq,  In  the  above  formula,  V represents  the  current  induced 
on  an  infinitely  long  cylinder  by  a unit  plane  wave: 


14ir  exp(-jkcos  0 *') 
V(0,z)  - - n 


H sin  0 W(kcos  0 ) * 

n n 


(3-7) 


in  which  W stands  for 


i(2) , 


W - - jtrj  (ka  sin  9 )H:r'  (ka  sin  9 ) 
0 w n 0 w n 


(3-8) 


,(2) 


where  and  Hq  are,  respectively,  the  zeroth-order  Bessel  function  of 

the  first  kind  and  the  Hankel  function  of  the  second  kind,  and  a is  the 

w 

radius  of  the  cylinder.  Another  universal  function  U is  found  in  the 
simple  formula;  it  represents  the  current  on  an  infinitely  long  center- 
fed  antenna  generated  by  a unit  voltage  impulse.  For  a thin-wire  antenna 
and  for  a sufficiently  large  kz,  U can  be  accurately  approximated  by 


u(z)  A il  . expCzjklL’-lJ 

UU;  n *n(2k  z'|)  - jir/2  - 2£n(ka  ) - y 

w 


(3-9) 


where  y ■ 0.57712...  is  thfe  Euler's  constant.  The  reflection  coefficient 
R in  Eq.  3-6  is  defined  by 

Jv0 


R(9,z)  - - ^^£n(kawsin  -|) 


+ Y + il  + * 


2 


•’] 


(3-10) 
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where 

VQ  - kz ' (1  - cos  0 ) , 

o n 

and  E^  is  the  exponential  integral.  Finally,  A and  C are, respectively, 
shorthand  notations  for 


A - R(9  »2h  ) - R(0  ,h  +z') 
n w n w 


(3-11) 


and 


_ _ R(t-0  ,2h  ) V (6.h  )R(ir,2h  ) U (2h  ) - R(0  ,2h  ) V (7r-0,h) 
C * n w n w w w n w n w 

1 - [R0r,2h  ) U (2h  )]2 
w w 


(3-12) 


(2) 

Again,  the  total  induced  current  is  obtained  by  superimposing  all  Ifl  • 

As  is  illustrated  in  Figure  3,  the  electric  field  of  the  incident  plane 

wave  is  in  the  0-direction  and  has  an  amplitude  f /sin  0 . Therefore, 

n n 

(2) 

I is  related  to  I by 
n 


I(z)  - l f (p,<|>)  -7 

"•  n c 


I^2)(z) 


sin  0 

n 


(3-13) 


where  P is  evaluated  at  the  location  of  the  wire. 

A careful  scrutiny  of  the  techniques  in  deriving  both  1^^  and  1^ 

reveals  that  the  two  methods  apply  only  if  the  incident  plane  wave  is 

2 

homogeneous,  i.e.,  cos  0^  < 1 or  real  incident  angle  ©n.  In  the  present 
problem,  however,  the  field  E inside  the  cavity  consists  of  both  homo- 
geneous and  inhomogeneous  components.  Hence,  both  methods  have  to  be 
extended  by  analytical  continuation  to  cover  the  case  of  an  incident 
inhomogeneous  plane  wave. 

Let  us  first  consider  1^^  which  is  obtained  by  solving  Eq.  3-4. 

Note  that  Eq.  3-4  is  developed  by  equating  the  axial  component  of  the 
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Figure  3.  A finite  length  cylinder  illuminated  by  an  incident  plane 
wave  at  an  oblique  angle. 


! 


electric  field  on  the  wire  and  that  of  the  incident  field.  The  tangential 

component  of  the  incident  electric  field  is  f exp(+1kcos  9 z)  which 

n n 

2 

remains  valid  even  in  the  case  when  cos  9 > 1,  Thus  no  modification 

n 

is  necessary  here,  although  allowances  in  the  actual  computer  programming 
must  be  made  to  accommodate  the  rapidly  oscillating  nature  of 
exp(+jkcos  9^z)  when  cos  9^  is  large. 

(2) 

Next,  we  have  to  extend  the  formula  of  I in  Eq.  3-6.  In  the  case 

n 

of  an  inhomogeneous  plane  wave,  sin  9 is  pure  imaginary,  and  may  be 

n 

calculated  from 


sin  9 * ±j  |./cos2  9-1  •,  cos2  9 > 1 

n I V n n 


(3-14) 


The  question  is  then:  "Which  sign  in  Eq.  3-14  should  be  used  in  the 
calculation?"  Consider  an  inhomogeneous  plane  wave  propagating  in  the 


x - z plane: 


-1k(x  sin  9 +z  cos  9 ) 
J n n 


; cos  9 > 1 

n 


(3-15) 


In  order  to  satisfy  the  radiation  condition,  the  field  must  decay  (instead 


of  grow)  exponentially  as  x + +°°.  This  imposes  a condition  on  sin  9^,  viz.. 


I (sin  9 ) < 0. 
m n 


(3-16) 


Thus,  the  lower  sign  (minus  sign)  in  Eq.  (3-14)  must  be  used  in  the 
(2} 

calculation  of  Iv  from  Eqs.  3-7,  3-8,  3-10,  and  3-13.  Furthermore, 
n 

W in  Eq.  3-8  becomes 

W - 2In (ka  Vcos2  9 - l)Kft(ka  i/cos2  9-1)  , (3-17) 

0 w’  n 0 w*  n 

where  1^  and  are,  respectively,  the  zeroth  order  modified  Bessel  functions, 


rids'lft* 


IV.  NUMERICAL  COMPUTATIONS 

As  it  stands,  the  numerical  computation  of  each  element  of  the  system 
of  linear  equations  in  Eq.  2-25  is  extremely  time-consuming.  The 
summation  of  the  ratio  of  Bessel  functions  often  presents  another  problem 
because  both  the  numerator  and  the  denominator  could  exceed  the  range  of 
the  computer  and  yet  the  quotient  is  still  not  small  enough  to  warrant 
the  termination  of  the  summation  process.  In  this  section,  a technique 
to  circumvent  the  above  difficulty  is  presented.  And,  at  the  same  time, 
it  significantly  improves  the  rate  of  convergence.  To  best  illustrate 
this  technique,  the  evaluation  of  is  discussed  in  detail. 

Besides  a multiplying  constant,  is  rewritten  below  with  the 


order  of  summation  and  integration  interchanged: 

2 


A 


(1) 

'uv 


r sin 
JQ  a' 


nA  00  gm,,cos  m<j)  H^2),(ya) 

Otd  r mV 0 id A 

l «■»  *•*  dot 


2Y  m*0  e (r2  - m2)  H^2)(Ya) 
my  m 


(4-1) 


Because  of  the  branch  point  a » k,  we  partition  the  above  integral  into 
two  parts: 


V(D 

UV 


k 

2 

sin  ad 

00 

v 

g COS 

6mv 

m<PQ 

• 0 

2 

a Y 

1 

m“0 

e (r2 

m u 

- m2) 

r 

sin2  ad 

00 

8mvCOS 

m4»0 

k 

a2t 

L 

m-0 

e (r2 

m u 

- m2) 

1(2)'( 

m 


<(2), 

m 


K'(ta) 
m 


da 


da 


(4-2) 


tli  / 2 2 

where  {K  } are  the  m order  modified  Bessel  functions  and  T ■ /a  - x . 
m 

Note  that  both  integrands  above  have  a non-integrable  singularity  in 
the  neighborhood  of  the  branch  point.  Hence,  both  integrals  must  be 
further  partitioned  such  that  the  integration  of  a small  region  around 
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the  branch  point  is  deliberately  Isolated.  Each  of  the  partial  integrals 
would  be  treated  individually  in  the  subsequent  paragraphs. 

The  first  integral  to  be  treated  is 


‘i-f 


k’S  IH  (T.)  da 

ay 


(4-3) 


where  6 is  an  arbitrarily  small  number  and  represents  the  summation  of 
Hankel  functions.  By  invoking  the  recurrence  relationship  between  the 
Hankel  function  and  its  derivative,  becomes 


(4-4) 


where  u • ya  and 


8mv“S"*0 


Gmv  /r2  2. 

m m 


m 


0 » 1 * 2 , • • 


(4-5) 


Since  the  real  and  imaginary  parts  of  ^ have  different  rates  of  convergence, 

they  are  treated  independently  as  follows:  First,  by  explicitly  writing 

as  J + jY  , the  real  part  of  Y reads 
in  in  in  ri 


Re{^H}  - -G( 


JqMJ^u)  + YQ(u)Y1(u) 

OV  _ 2 / \ j.  V2/  \ 

J0(u)  + Yq(u) 


+ l « 


m“l 


mv 


J (u) J .(u)  + Y (u)Y  . (u) 
m m— 1 m m— l 


J2(u)  + Y2(u) 
m in 


(4-6) 


Note  that  for  a given  u,  the  large-order  approximations  of  each  trrm  of 
the  3um  can  be  extracted,  which  are 


lim 

nr*=° 


J J . + YY  , 
m in- 1 m m-i  m 


2 2 

J + Y 
m m 


m 

u 


(4-7) 
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and 


sin  2m4>r 


lim  Gnv  * " cos  v*0 


(4-8) 


m 


By  adding  and  subtracting  the  large-order  approximations  from  each  term 
of  the  sum,  Eq.  4-6  becomes 


Re{£H}  - - G( 


Ov 


Vl  * Vl 
2 2 
Jo  + Yo 


jj  ,+n  . 

p nr"  1 in  ni~  1 tn 

J2  + Y2  U 

m m 


cos  v<(>Q  sin 


$ 


cos  v<f>-  « sin  2m$. 

+ l 5— 2 

u L.,  2 

tn-1  m 


The  second  sum  can  be  evaluated  analytically  and  the  result  is 


(4-9) 


® sin  2mt>  rzq> 

l ^ - ° ln(>2  sin  f>  dt 
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The  integral  above  is  a thoroughly  studied  special  function,  known  as 

Clausen's  integral  [14],  whose  value  can  be  easily  determined.  The 

remaining  sum  in  Eq.  4-9  has  to  be  evaluated  numerically.  However,  as 

-2 

compared  with  the  sum  in  Eq.  4-6,  which  converges  at  a rate  of  m , the 

-4 

modified  sum  converges  at  a much  faster  rate  of  m . Furthermore,  in 

evaluating  the  m*"*1  order  Bessel  function  B (B  ■ J or  Y ) , we  can 

m m m m 


apply  the  following  recurrence  formula: 


2m 


(4-11) 


B (u)  « — B (u)  - B ,(u)  . 

m+1  u m m-1 

However,  as  is  well-known,  we  should  use  this  recurrence  formula  with 

extreme  caution  in  computing  J to  avoid  the  so-called  "propagation  of 

P 

error"  when  2m/u  >1.  In  the  actual  computation,  the  total  sum  is  broken 

down  into  partial  sums  of  10,  e.g. , from  n to  n + 9.  We  first  evaluate 

J and  J then,  we  apply  the  recurrence  formula  in  both  the  forward 
n+4  n+5 


and  backward  directions  to  obtain  the  rest  of  J . Moreover,  we  terminate 

m 
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the  sunznation  by  comparing  the  magnitude  of  the  partial  sum  with  the  total 
sum.  In  doing  so,  we  avoid  the  danger  of  terminating  the  summation  pre- 
maturely in  the  case  2n<j>g  is  a multiple  of  if  (.referring  to  Eq.  4-8). 

As  a final  remark,  when  ka  » 1,  L can  be  approximated  by 

w n 

!„<»>  i ^ Jg(u)  ; „ > 100  . (4-12) 

On  the  other  hand,  no  modification  is  necessary  to  sum  Im{£H}  because 
its  large-order  approximation  is  zero. 

The  second  integral  of  concern  is 
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where  £ represents  the  summation  of  the  modified  Bessel  functions  {K^} 
With  the  derivative  of  K written  in  terms  of  and  K^,  reads 
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in  which  v - Ta.  Since  the  large-order  approximation  of  each  term  of  the 
sum  is 
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we  evaluate  as 
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The  above  formula  is  similar  to  Eq.  4-9;  therefore,  it  is  evaluated  by 
similar  techniques.  Furthermore,  we  extract  the  large  argument  approximation 
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of  £ , which  is 


L(v)  : L(150)  ; V > 150  . 


(4-17) 
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Substituting  the  results  of  Eqs.  4-16  and  4-17  into  Eq.  4-13,  we  arrive  at 
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By  two  successive  integrations  by  parts,  the  last  integral  above  is 
transformed  into 
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where  C ^ is  the  cosine  integral  [15]. 
The  third  integral  of  concern  is 
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First,  we  examine  the  small- argument  behaviors  of  both  and  VR;  they  are 
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It  can  be  easily  shown  that  the  real  part  of  I3  tends  to  zero,  and  its 
imaginary  part  can  be  approximated  by 


ttG„  s in  kd 
1 ~ 
m 3 2 ax  £n(a/2kj) 


(4-22) 


Since  the  forms  of  and  A^^  are  similar,  analogous  techniques 

are  employed  to  evaluate  A^  ; the  result  is  (besides  a multiplying  constant) 
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where  N » Integer  (kh/ir)  and 
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The  last  sum  in  Eq.  4—23  can  be  determined  as  follows: 
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where  f is  the  Clausen's  integral  stated  in  Eq.  4-10. 
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(4-25) 
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V.  NUMERICAL  RESULTS 


In  determining  the  field  inside  the  cavity  with  the  wire  absent, 
the  series  representation  of  the  aperture  field  in  Eq.  2-3  is  truncated 
at  n ■ N ( the  series  contains  N + 1 terms) . We  must  first  establish  the 
convergence  of  the  aperture  field  with  respect  to  N.  The  aperture  fields 
as  calculated with  N ■ 0,  1,  and  2 are  shown  in  Figure  4.  We  notice  that 
the  aperture  fields  as  computed  by  N ■ 1 and  that  by  N ■ 2 agree  reasonably 
at  the  main  lobe  but  not  so  well  at  the  side  lobe.  For  many  practical 
cases,  the  three-term  expansion  N = 2 is  generally  sufficiently  accurate 
for  computing  the  aperture  field. 

An  indication  of  the  accuracy  of  the  field  E inside  the  cavity  is 
how  well  does  E satisfy  the  boundary  condition  on  the  cavity  wall  (including 
the  aperture).  In  Figure  5,  E as  computed  from  Eq.  2-26  is  plotted  as  a 
function  of  4> . In  the  aperture  defined  by  |$(  £ 57.29°,  the  calculated 
Ez  agrees  extremely  well  with  the  two-term  expansion  of  the  aperature  field 
calculated  from  Eq.  2-3.  It  drops  to  less  than  0.01%  of  the  aperture 
field  on  the  wall  where  E^  should  be  ideally  zero.  Also  shown  in  Figure  5 
is  E^  at  points  just  behind  the  aperture,  p/a  * 0.995.  The  variations 
of  E^  with  respect  to  z at  p/a  ■ 0.0,  0.5,  and  0.8  are  sketched  in  Figure  6. 
As  a function  of  decreasing  (p/a),  Ez  decreases  rapidly  from  its  value 
in  the  aperture,  while  it  increases  at  an  even  faster  rate  from  zero  on 
the  cavity  wall  toward  the  center  of  the  cavity.  These  features  are 
illustrated  in  Figure  7 where  Ez  is  plotted  as  a function  of  p at  z * 0.0 
and  z/h  ■ 0.3. 

Part  B of  the  problem  is  to  compute  the  current  induced  on  the  wire 
inside  the  cavity.  As  mentioned  in  previous  sections,  the  simple  formula 
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Figure  4.  E_  in  the  aperture  as  a function  of  <t>  with  N as 

a^parameter.  The  input  data  ares  a **  0.3X,  c = 0.3X, 
d - 0.015X,  h - 0.6X. 
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Figure  5.  E inside  the  cavity  as  a function  of  <J>  with  (p/a) 
as  a parameter.  The  input  data  are:  a * 0.3X, 


E.(z)/E 


z/h 


Figure  6.  Ez  inside  the  cavity  as  a function  of  z with  (p/a) 
as  a parameter.  The  input  data  ares  a * 0.3X, 
c ■ 0.3X,  d ■ 0.015X,  h ■ 0.6X,  * 0. 
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Figure  7.  Ez  inside  the  cavity  as  a function  of  p with  (z/h) 
as  a parameter.  The  input  data  are:  a » 0.3X, 
c - 0.3X,  d - 0.015X,  h - 0.6X. 


of  che  length  of  the  slot  2c.  The  induced  current  at  the  center  of  the 
wire  does  increase  with  c when  c is  relatively  small  as  compared  with 
the  radius  of  the  cavity  a.  The  center  current  acquires  its  peak  value 
when  c/a  is  about  0.9  and  begins  to  decrease  thereafter.  These  charac- 
teristics are  demonstrated  in  Figure  11. 

In  Figure  12 , we  show  che  center  current  on  the  wire  as  a function 
of  a,  0.25  £ a/A  £ 0.8,  with  the  other  dimensions  fixed.  We  observe  that 
there  are  three  peaks.  The  first  of  these  peaks  corresponds  to  the 
resonance  at  which  J^y^a)  * 0.  The  second  and  third  peaks  correspond 
respectively  to  J^Cy^a)  * 0 and  Jg(y^a)  » 0.  It  should  be  noted  chat  the 
resonances  would  occur  more  often  beyond  the  range  of  Figure  12,  and 
numerical  solutions  for  large  a would  not  be  reliable. 

The  center  current  on  a wire  as  a function  of  the  length  of  the 
cavity  2h  with  other  dimensions  fixed  is  plotted  in  Figure  13.  Two 
peaks  are  observed  in  Figure  13.  The  first  peak  is  related  to  the 
resonance  corresponding  to  J^y^a)  * 0,  but  we  are  unable  to  analyze 
' the  nature  of  the  smaller  peak  at  h * 1.02A. 

Finally,  we  examine  the  variation  of  the  induced  current  with  respect 
to  the  frequency  of  the  incident  plane  wave.  The  induced  currents  at 
representative  frequencies  are  sketched  in  Figure  14,  which  clearly  reflects 
the  alteration  of  the  electrical  length  of  the  wire. 
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Figure  11.  Induced  current  at  the  center  of  the  wire  inside  the  cavity 


as  a function  of  slot  length  2c.  The  input  data  are: 
a * 1,  d ■ 0.015,  h - 2.2,  freq.  - 0.3. 
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a/X 

Figure  12.  Induced  current  at  the  center  of  the  wire  inside  the 
cavity  as  a function  of  cylinder  radius  a.  The  input 
data  are:  c ■ 0.3A,  d » 0.015A,  h » 0.6A,  h ■ 0.2A, 

U 

p/a  - 0.1. 
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Figure  13.  Induced  current  at  the  center  of  the  wire  inside  the  cavity  as  a function  of  cavity 
length  2h.  The  input  data  are:  a = 0.5A,  c = 0.3A,  d = 0.015A,  h = 0.2A, 
a = 0.001A,  p/a  = 0.1.  w 
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Figure  14.  Current  distribution  induced  on  a wire  inside  the 

cavity  with  frequency  as  a parameter.  The  input 

data  are:  a ■ 0.35m,  c » 0.25m,  d « 0.005m, 

h * 0.82m,  a ■ 0.001m. 
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ABSTRACT 


The  problem  of  radiation  from  sources  in  the  presence  of 
smooth,  convex,  impenetrable  objects  is  considered,  and  a brief 
survey  of  various  high  frequency  techniqes  is  presented.  A 
generalization  of  the  geometrical  theory  of  diffraction,  and  two 
new  techniques  based  on  the  spectral  domain  approach  and  an 
asymptotic  evaluation  of  the  radiation  integral  for  the  surface 
current,  also  are  discussed.  Some  numerical  results  derived  from 
the  spectral  domain  formulas  are  presented  and  a comparison  with 
available  theoretical  and  experimental  data  is  included. 
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presumably  attributable  to  experimental  error  36 
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1.  Introduction 


The  problems  of  radiation  from  sources  in  the  presence  of  impenetrable 
smooth  convex  objects  and  the  diffraction  of  a plane  wave  by  such  objects  are 
of  great  practical  interest  in  the  design  of  antennas  on  structures,  e.g., 
conformal  arrays.  Unfortunately,  the  exact  analytical  solutions  to  these 
problems,  based  on  the  methods  of  "separation  of  variables"  or  "function- 
theoretic"  procedures  (Wiener-Hopf  technique,  residue  calculus,  etc.),  exist 
only  for  a very  limited  number  of  scattering  geometries.  Furthermore, 
the  exact  solutions  are  typically  highly  complex  in  nature;  hence,  the  process 
of  extracting  numerical  results  from  them  can  be  very  time-consuming  and  is 
by  no  means  trivial.  This  situation  has  motivated  many  researchers  to  explore 
approaches  to  the  problems  of  radiation  and  scattering  from  smooth  convex 
. structures. 

In  the  low  and  resonant  frequency  ranges,  several  reliable  numerical 
procedures,  e.g.,  the  moment  method,  are  available  for  solving  the  radiation 
and  scattering  problems.  However,  in  the  high  frequency  domain,  numerical 
techniques  based  on  matrix  methods  become  unwieldly  if  noc  impractical,  prompting 
one  to  employ  asymptotic  techniques  suitable  for  large  k(*2i/\),  where  \ is 
the  wavelength  of  the  illuminating  wave. 

In  this  work,  we  begin  by  presenting,  in  Sec.  2 , a survey  of  various 
high  frequency  asymptotic  techniques  for  the  problem  stated  above.  The  survey 
will  be  necessarily  brief,  and  will  cover  only  the  highlights  of  a number  of 
important  approaches  to  the  problem  at  hand,  viz.,  Fock's  theory,  the  geometrical 
theory  of  diffraction  (GTD) , and  the  direct  integral  equation  approach.  The 
reader  interested  in  further  details  may  choose  to  consult  the  works  of  Bowman, 
at  al.  [1],  Uslenghi  [2],  and  Xouyoumjian  [3]. 
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In  Sec.  3,  we  consider  the  generalization  of  GTD  and  present  3ome  new 
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approaches  to  the  curved  surface  radiation  and  scattering  problems.  Some 
numerical  results  based  on  one  of  these  new  approaches  are  presented  in  Sec.  3 
and  a comparison  with  other  available  methods  are  Included. 

2.  Survey  of  Available  High-Freauencv  Asymptotic  Techniques 
2.1  Watson  Transformation 

One  of  the  first  successful  attempts  to  derive  an  asymptotic  expansion 
for  the  far-field  generated  by  a point  source  located  in  the  proximity  of  a 
conducting  surface  was  made  by  G.  N.  Watson  in  1913  [4].  His  method, 
essentially,  consisted  of  two  steps:  1)  transforming  the  original  infinite 
series  solution  into  a contour  integral  (by  Cauchy's  residue  theorem);  2) 
deforming  the  contour  of  integration  30  as  to  capture  a set  of  complex  poles 
of  the  integrand.  The  original  integral  is  then  expressed  in  terms  of  an 
infinite  series  which  converges  very  rapidly,  provided  the  observation  point 
is  in  the  shadow  region.  The  first  few  terms  of  this  series  were  later  inter- 
preted as  "creeping  waves."  The  method  was  first  applied  to  a sphere  and 
circular  cylinder,  and  later  to  some  other  geometries  as  well.  The  mathematical 
rigor  of  the  method  was  the  subject  of  further  investigations  by  other  researchers 
([5],  [6]  and  [7]).  Although  Watson  transformation  can  only  be  applied  to  a 
few  simple  goemetries,  e.g.,  the  sphere,  cylinder,  cone,  spheroid,  it  is 

still  regarded  as  one  of  the  cornerstones  of  the  more  general  high  frequency 
techniques  because  of  its  mathematical  rigor.  Watson  transformation  is 
especially  powerful  in  the  shadow  region  of  the  geometric  optics  field.  In 
the  lit  region,  the  above-mentioned  contour  integral  is  evaluated  using  the 
"stationary  phase"  method  and  yields  the  reflected  field  from  the  surface. 

In  this  region,  the  most  significant  contribution  to  the  total  scattered  field 
typically  comes  from  the  surface  currant  induced  on  the  smooth  convex  part  of 
the  object;  the  so-called  "Physical  Optics"  approximation  can  be  applied 
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([3],  [9],  and  [10])  to  derive  the  reflected  field.  The  Physical  Optics 
method  is  based  upon  approximating  the  induced  surface  currenc  in  the  lit 
region  of  the  object  by  the  current  that  would  be  induced  on  the  local  tangent 
plane,  and  by  assuming  that  the  surface  currenc  is  zero  in  the  shadow  region. 

The  far  field  is  constructed  by  substituting  the  above  estimate  for  the 
induced  surface  currenc  in  the  integral  representation  of  the  scattered  field, 
and  evaluating  the  same  in  an  asymptotic  sense.  The  dominant  term  of  the 
asymptotic  expansion  of  this  integral  can  be  shown  to  be  identical  to  the  first  term 
of  the  Luneberg-Kline  expansion  of  the  geometrical  optics  far  field  ([11]  and 
[12]).  However,  the  higher-order  terms  derived  from  the  physical  optics 
approach  do  not  provide  us  with  the  correct  result  in.  the  shadow  or  transition 
regions  where  the  diffracted  field  contributes  the  most. 

In  the  next  section,  we  discuss  Fock's  theory,  which  can  fill  the  gap 
between  the  Physical  Optics  in  the  lit  region  and  the  "creeping  wave" 
representation  in  the  shadow  region. 

2. 2 Fock's  Theory 

The  region  between  the  lit  and  the  shadow  part  on  a surface  is  called 

"penumbra  region."  The  angular  width  of  this  region  is  approximately  given 
, 1/3 

by  (aT“/t)  where  A is  the  wavelength  of  the  illumination  and  r is  the 
o o 

radius  of  curvature  of  the  surface  of  the  object  in  this  region  in  the  incident 
plane  (Fig.  1).  Fock's  theory  invokes  the  principle  of  Icoal  character  of 
the  field  in  the  penumbra  region  [13]  and  is  based  on  the  conjecture  that  all 
bodies  with  a smoothly  varying  curvature  have  the  same  current  distribution 
in  the  penumbra  region,  provided  that  the  curvature  and  the  incident  wave  are 
the  same  near  the  point  under  consideration.  This  principle  allows  one  to 
locally  replace  the  surface  of  the  object  by  a portion  of  a paraboloid  of 
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Figure  1:  Section  of  cha  body  ia  “he  plane  of  incidence. 
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?•*  rure  2:  Plane  wave  incident  uoon  a smooth  convex  body. 
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revolution.  A unique  feature  of  the  expressions  for  Fock  currents  is  that 
they  provide  a convenient  transformation  of  the  geometric  optics  currents 
in  the  lit  region  into  the  creeping  wave  currents  in  the  shadow  region. 

Fock  himself  deduced  Che  pertinent  formulas  for  the  surface  currents  by 
treating  a convex  body  problem  [14]  described  below. 

Consider  a convex  body  and  a plane  wave  incident  in  the  direction  of  the 
x-axis.  If  the  equation  of  the  surface  is 


f(x,y,z)  - 0 


(1) 


then  the  curve  representing  the  boundary  of  geometrical  shadow  is  given  by 


f(r) 


0 * If 


(2) 


Consider  a point  0 on  the  boundary  of  a shadow  region  where  we  set  up 
a rectangular  coordinate  system  as  shown  in  Fig.  2 (z:  normal  to  the  surface, 
x:  in  the  direction  of  propagation,  and  y is  the  tangent  to  the  boundary  of 
shadow).  In  the  vicinity  of  this  point,  the  surface  of  the  body  could  be 
locally  replaced  by  a paraboloid  of  revolution  which  is  expressed  by  the 
equation. 


2 2 
z + 1/2  (ax“  + 2bxy  + cy  ) 


(3) 


Each  of  the  field  components  satisfies  the  Helmholtz  equation 


2 2 

(7  + k ) H - 0 


(4) 
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The  face  chac  the  incident  wave  travels  along  the  x-axis,  suggests  that  u. 
be  written  in  the  fonn 


* . ie-Jkx 


where  an  exp(j^e)  time  dependence  has  been  assumed.  Substituting  (5)  in  (4) 


gives 


72?  - 2Jk  ^ - 0 


At  this  point,  two  basic  assumptions  are  introduced  in  Fock’s  theory, 


i)  V's  are  relatively  slowly  varying  function  of  coordinates 

ii)  ? varies  more  rapidly  in  the  z-direction  than  in  x and  y, 

i.e.. 


XV  k - k - k~ 

b m *> . hm  > b - 0 £ 

3z  m 3x  m 3y  m 


Based  upon  (7),  we  can  write  (6)  as 


3 v 3 -v 

-2— i-  - 21k  — » 0 
3zz  J 3x 


2 2 
and  consequently  m'  ■ m"  (m  i3  very  large),  where  the  terms  of  order  1/m 

have  been  omitted. 

Inserting  these  estimates  and  assumptions  into  the  Maxwell's  equation, 
we  can  find  seme  simple  expressions  for  all  the  field  components  in  terms  of 


H and  3 . If  we  write  3 as 

y z y 


(9) 


H - H°  «’jkX  J 

7 7 


where  H°  is  che  magnitude  of  the  incident  wave  at  infinity,  then  '?  must 
satisfy 


C1Q) 


with  boundary  condition 


3Y 

3z 


_1 

-jk  (ax  + by  ) 


0 


(ID 


on  the  surface  of  the  body.  Eqn.  (11)  is  the  simplified  version  of  the 
Leontovich  boundary  condition  where 


n ■ e - j 


4iro 

ek 


The  final  solution  for  H on  the  surface  of  the  body  which  satisfies 

y 

the  boundary  condition  and  the  condition  at  infinity,  may  be  written  in  the 
form 

H - Hex  G(C,q)  (12) 

7 7 


where  H6*  » external  field 


G(5,q)  - e"(J/3)  5 V1(5,q) 


V^(^,q)  ■ Fock  function  defined  in  the  Appendix  A. 

* m(ax  + by)  ■ reduced  distance  from  the  shadow  boundary  * <./d. 

2 1/3 

d • the  width  of  penumbra  region  ■ (2rQ  /k) 

i • distance  between  che  observation  point  and  the  shadow  boundary  along  the 
incident  ray  (Fig.  3). 
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Figure  4:  Coordinates  of  observation  point  in  terns  of 


q * -jm//n-  -(j/n)  \rh  (■  0 for  conducting  body) 


The  ocher  tangential  component  of  the  magnetic  field  11^  on  the  surface 
of  the  body  can  be  obtained  in  a similar  manner 


H - H**  [-  -1  e-(j/3)C  /3  f(£)  ] 

x z m 


(13) 


where  f(5)  is  another  Fock  function  defined  in  Appendix  A.  Fock's 
formulas  not  only  give  the  surface  value  of  the  field,  but  also  can  be 
utilized  to  find  the  field  in  the  proximity  of  the  object.  For  a plane  wave 
incidence,  the  first  order,  i.e.,  0(l/m)  terms  for  the  scattered  field  within 
a certain  layer  around  the  object,  can  be  written  as 


Hx«°,  Hy  - e"jkX  V («,«,  ^ • H°  e‘jkx  *}(r,n) 


E -(j/m)H°  3Y/3C,  E - H , E - -H 

X y v 2 2 v 


(14) 


where 


2.2  2 

C ■ 2am  [z  +(l/2)(ax“  + 2bxy  + cy  )J  * reduced  heighc  from  the  surface 


of  the  body  (Fig.  4) 


w^(t)-qw1(t) 


- v (t) 


(15) 


The  path  of  integration  for  $ and  ? is  shown  in  Fig.  5 


dc 
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Complex  t*  plane 


*c«-  ■■ 


Figure  5:  Path  of  integration  for  (15)  in  the  complex  t-plane 


Figure  5:  Comparison  between  the  various  definitions  of 

parameter  % for  the  case  of  a circular  cylinder. 
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Fock's  expressions  for  Che  field  components  in  che  penumbra  region 
(C  - 0)  can  be  extended  to  Che  shadow  region,  by  introducing  some  modifications 
in  the  definition  of  parameter  £.  Goodrich  [15]  has  generalized  the 

argument  used  by  Fock  in  the  penumbra  region  to  anywhere  in  the  shadow 
region  by  introducing  a new  set  of  variables,  £ and  C,for  che  incremental 
distances  along  the  path  leading  into  the  shadow  region.  In  this  generalization, 
the  parameter  • as  defined  in  (12)  is  replaced  by 


ds 

R(s) 


(16) 


where  s is  the  arc  length  along  the  geodesics  which  originate  from  che  shadow 
boundary  and  go  into  the  shadow  region  along  the  surface,  and  R(s)  is  che 
radius  of  curvature  of  the  surface  along  che  geodesics.  For  the  case  of  a 
circular  cylinder  of  radius  a (Fig.  6),  the  expression  C simplifies  to 


5 - (ka/2)1/3  9 - s/d  (17) 

Fock  also  treated  the  case  where  the  point  source  was  very  close  to 
the  surface  of  the  body.  He  analyzed  the  radiation  of  electric  dipoles  near  a 
spherical  model  of  the  earth  [16]  and  derived  the  formulas  for  the  scattered 
fields  in  terms  of  functions  (attenuation  functions)  similar  to  p and 
which  are  valid  both  in  the  shadow  and  transition  regions  [17].  Fock's 
assumptions  were  later  proven  in  a more  systematic  and  mathematically  rigorous 
manner  by  Cullen  [18]  and  Hong  [19]  by  using  a direct  integral  equation 
approach.  This  method  is  described  in  the  next  section. 
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2. 3 Direct  Integral  Equation  Approaches 

This  method,  which  is  closely  related  to  Fock's  theory,  can  be 
illustrated  by  analyzing  the  diffraction  of  a plane  electromagnetic  wave 
by  an  arbitrary  conducting  body  (large  compared  with  \) . Cullen  [13]  obtained 
a first-order  asymptotic  solution  to  the  integral  equation  for  the  induced 
surface  current 


J(r)  - 2n(r)  * Hinc(?)  - (l/2ir)n(r) 


f / ds'  tf(r’)  , 

S R 


where  n(r) is  the  outward  unit  normal  to  the  surface  at  r,  Hinc(r)  is 
the  incident  magnetic  field  on  the  surface  (S)  of  the  body, and  R - r - r' 

(r ' is  a variable  point  on  the  surface) . 

Fock  used  this  integral  equation  to  deduce  the  important  principle 
of  local  character  of  the  field  in  the  penumbra  region.  Cullen  derived  a 
first-order  asymptotic  solution  to  (18)  which  agreed  with  Fock's  results 
given  in  (12)  and  (13).  Cullen's  method  consists  of  transferring  the  two- 
dimensional  integral  equation  (13),  in  the  penumbra  region,  to  a one-dimensional, 
Volterra-type  equation.  This  is  accomplished  by  applying  the  stationary 
phase  technique  to  the  original  integral  while  integrating  with  respect  to 
one  of  the  variables.  The  resulting  one-dimensional  Volterra  equation  is  then 
solved  in  Cullen's  method  by  the  Fourier  transform  technique.  A similar 
procedure  was  used  by  Hong  [19]  to  analyze,  asymptotically,  the  diffraction 
of  electromagnetic  and  acoustic  plane  waves  by  smooth  convex  bodies.  We  will 
now  proceed  to  explain  Hong's  method  in  a little  more  detail  by  referring 
back,  once  again,  to  the  integral  aquation  (IS).  The  surface  is  parametrized  by 
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che  geodesic  coordinate  system  (o,v)  such  chac  che  shadow  boundary  for  the 
incident  plane  wave  traveling  along  che  tangent  o(0,v)  to  the  v ■ 0 curve  is  che 
o « 0 curve.  The  quantities  o(o,v),  b(o,v)  and  n(o,v)  form  a right-hand  local 

A A ^ 

orthonormal  basis  (n  ■ a * b)  (Fig-  ') • 

Since  the  incident  field  has  a phase  factor  e jk°(0,0)  r(u,0)^ 
we  write  the  surface  current  in  the  form 


J(r)  - [ (r)  a(r)  + I^Cr)  Mr)]  e 


-jka 


(19) 


where  a is  the  arc  length  along  che  geodesic.  Substituting  (19)  back  into 
(18)  and  restricting  che  resulting  equation  to  che  points  on  the  geodesic  v-0, 
we  obtain  two  coupled,  two-dimens ional  integral  equations  for  l^(a,0)  and 
I, (s,0).  It  can  be  shown  that  these  integrals  have  saddle  points  at  v»0  (for 

D 

the  v- integration) . Applying  the  "steepest  descent  path"  method  to  v- integration, 


and  keeping  che  terms  up  to  the  order  1/M“,  where  Ji- 


ao^ (a , 0) ) 1/3 , we 


*0’  0 -'-cj 

obtain  the  following  decoupled  one-dimensional, Volterra-type  integral  equations 
for  I.  U,0)  and  I (5,0) 

D y 


Iff(5,0)  - 2 I*nc(?,0)  - / dr  Ia(5,0)  Ka(--t)  +0(MqJ) 


-3. 


(20) 


Ib(C,0)  - 2 Ibnc  (5,0)  - / dt  Ib(5,0)  ^(C-t)  + 0(M~J) 


,-3, 


3„(J»V)  is  the  radius  of  curvature  of  che  surface  along  geodesics  (v  * constant 
curves)  ac  point  (o,v). 

Solving  (20)  by  Fourier  transforms,  we  obtain  the  expression  for  che 
induced  currents  in  che  penumbra  and  shadow  regions, and  the  first-order  solutions 
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are  found  co  be  the  same  as  those  of  Fock  and  GTD  [20].  One  of  the  important 

conclusions  drawn  from  Hong's  solution  is  that  the  leading  term  in  the 

asymptotic  expansion,  which  is  the  same  for  the  acoustic  and  electromagnetic 

problems,  is  independent  of  curvature  in  the  direction  transverse  to  the 

geodesic,  provided  the  divergence  factor  is  suppressed.  However,  we  should 

bear  in  mind  that  Hong's  method  was  designed  for  the  case  of  axial  Incidence 

on  symmetric  objects, and  in  this  case,  the  geodesics  are  torsionless.  The 

above  conclusion  does  not  seem  to  be  valid  in  the  cases  where  the  rays  have 

nonzero  torsion  ([21],  [22]).  In  Hong's  expressions  for  the  surface  current, 

the  transverse  curvature  has  only  a second-order  effect.  It  was  also  shown 

-2/3 

that  up  to  the  terms  of  order  (ko^)  in  the  asymptotic  expansion,  the 
tangential  and  binormal  components  of  the  creeping  waves  are  not  coupled. 

Both  Fock's  theory  and  the  "direct  integral  equation  approach"  give  the 
induced  surface  current,  or  the  scattered  field  in  the  neighborhood  of  the 
surface  of  the  scatterer,  due  to  an  incident  plane  wave.  These  expressions  can 
also  be  used  to  derive  the  radiated  field  via  the  use  of  the  reciprocity  theorem. 

The  methods  which  have  been  discussed  thus  far  are  mathematically 
rigorous.  However,  they  are  limited  in  the  scope  of  their  application  to 
geometries  satisfying  some  special  smoothness  and  symmetry  criteria.  "Geometrical 
theory  of  diffraction"  (GTD),  which  we  discuss  in  the  next  section,  has  a 
broader  scope,  although  it  does  lack  the  mathematical  rigor  of  approaches 
described  until  now. 

2.4  Geometrical  Theory  of  Diffraction  (GTD) 

Geometrical  theory  of  diffraction  (GTD),  developed  by  J.  3.  Keller 
([20],  [23],  [24],  [25],  and  [26]),  is  a generalization  of  geometrical  optics. 
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It  is  based  upon  the  assumption  that  fields  propagate  along  rays.  Keller's 
major  contribution  was  to  introduce  the  new  kinds  of  rays  called  the  "diffracted 
rays"  which  together  with  the  geometrical  optics  rays, constitute  the  total 
field.  In  our  problem,  viz.,  source  radiating  in  che  proximity  of  the  smooch 
object,  che  diffracted  rays  travel  along  the  curves  on  che  surface  of  the 
scatterer.  By  applying  Fermat's  principle  to  these  surface  rays,  we  conclude 
that  the  above-mentioned  curves  should  be  geodesics  on  the  surface  of  the  body. 

In  che  GTD  procedure,  one  assigns  a value  to  che  field  along  each  ray  of  these 
surfaces.  The  total  field  at  any  point  in  che  space  is  the  sum  of  che  fields 
due  to  various  rays  (incident,  reflected  and  diffracted)  passing  through  that 
point.  An  important  advantage  of  che  GTD  approach  is  that  it  can  be  applied 
to  both  scalar  (acoustic)  and  vector  (electromagnetic)  problems  and  to  smooth 
convex  objects  of  an  arbitrary  shape. 

Consider  the  problem  of  determining  the  radiated  field  of  a scalar  point 
source  located  on  the  surface  of  a smooth  convex  opaque  body.  If  the  observation 
point  is  in  the  shadow  region,  the  ray  paths  originating  at  Q and  reaching  P 
(observation  point)  are  comprised  of  two  sections.  One  of  these  sections  follows  the 
straight  line  path  P^P,  while  che  other  travels  along  a geodesic  on  the  surface 
(Fig.  3).  Let  us  consider  the  propagation  of  the  field  along  each  section  separately. 

a)  Rays  in  free  space:  3ehavior  of  the  fields  along  these  rays  can  be 
determined  by  obtaining  a high-frequency  asymptotic  solution  to  Maxwell's 
equacion  ina source-free  homogeneous  isotropic  medium.  We  begin  with  the 
Luneberg-Kline  asymptotic  expansion  of  the  electric  field  ([11]  and  [12]): 


and  inserc  it  into  the  Maxwell's  equations.  This  results  in  the  following 


equations  governing  the  propagation  of  electromagnetic  fields  along  Che  rays. 


[7S(r)]  - 1 


(Eikonal  equation) 


2(7S*7)e  + (7zS)e  ■ -74  e , (Transport  equation) 

m m o-i 


7S ■ e ^ -7- e , (Gauss's  Law) 
m m-1 


e ^ * 0,  m ■ 0,  1,  2, 


The  2eroth-order  solution  to  Che  above  system  of  aquations, which  turns  out 


to  be  in  agreement  with  what  one  would  obtain  by  geometric  optics,  may  be 


written  as 


E(cr)  - EQ(0)e 


-Jlc0S(0) 


(pj+o)  (p  ) 


where  3 is  the  distance  traveled  along  the  ray  from  the  reference  point 


0(c«0)  on  the  ray  path,  p^  and  o^  are  che  principal  radii  of  curvature  of 


the  wavefront  at  3-0.  It  is  apparent  that  the  expression  fails  when  3«-o^  or 
a * -p^,  t.e.,  at  the  caustic  lines  (Fig.  9).  In  the  cases  where  it  is 


convenient  to  choose  the  point  of  diffraction  on  the  surface  of  a body  as 


the  reference  point  0,  the  formula  (2.5)  should  be  modified  as  follows 


'0  * \/  3(p+a) 


Ik 
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In  these  cases,  the  point  of  diffraction  itself  is  a caustic,  and  o is  the 
distance  betveen  this  point  and  the  second  caustic. 

b)  Surface  Rays:  These  rays  follow  the  surface  S along  the  geodesics 
into  the  shadow  region,  and  shed  off  energy  tangentially  as  they  propagate. 

In  order  to  study  the  behavior  of  the  field  along  these  rays,  we  introduce  a 
special  ray-fixed  coordinate  system, a, n,b. 

a:  Unit  vector  tangent  to  the  ray;  n:  outward  unit  normal  to  the 
surface;  and  b * t * n or  binormal  direction;  a vector  field  can  be  decomposed 
into  its  components  along  these  unit  vectors  as 

E - E^a  -t-  E^n  + Ebb  (27) 

At  this  point,  several  important  assumptions  are  introduced  in  the  GTD 
approach  [20]: 

i)  E and  H are  orthogonal  to  each  other  and  to  the  ray. 

ii)  Variation  of  the  phase  of  the  field  along  the  ray  is 
the  same  for  both  fields. 

iii)  E^  and  E^  propagate  independently,  and  E^  ■ 0. 

iv)  satisfies  the  scalar  wave  equation  (V2+k2)  u«0  with  the 

boundarv  condition  u«0  on  the  surface  S,  while  E satisfies 

n 

the  same  equation  with  the  boundary  condition  * 0. 

on 

The  next  step  in  the  GTD  approach  is  to  conj ecture, on  the  basis  of  the 
solution  to  some  canonical  problems,  that  the  surface  field  propagating  along 
each  ray  is  comprised  of  an  infinite  set  of  "modes."  Along  a ray-fixed  path 
GTD  assigns  a complex  value  to  each  component  of  the  field  associated  with 
the  individual  modes.  The  propagation  of  these  modal  field  is  described  by 
the  equation 
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j (*0-ico) 


a(o-)  - A(c)  e 


C28) 


when  a is  che  distance  between  an  arbitrary  point  along  the  ray  and  the  source 
Q and  ^ is  the  phase  of  Che  field  at  che  source  point.  Next,  invoking  the 
principle  of  conservation  of  energy  between  two  adjacent  ray3,  and  using  the  fact 
chat  che  surface  rays  shed  energy  off  tangentially,  we  can  arrive  at  the 
following  expression  for  a(o) 


exp [ - j ka 


a 


J * 


0 


(o' )do' ] 


(29) 


where  a (a)  is  the  "attenuation  constant,"  K is  proportional  to  the  strength  of 

1/2 

che  source,  and  df^,  and  c are  shown  in  Figure  10.  The  quantity  [d*^/  (pd'?^)  1 
indicates  the  spreading  of  *e  surface  ray  tube"  as  it  travels  along  the  surface. 
Equations  (26)  and  (29)  describe  the  laws  of  propagation  for  the  rays  which 
originate  from  the  source  point  Q,  are  diffracted  at  ?^,  and  reach  the 
observation  point  P.  To  complete  the  solution,  we  need  to  determine  the  actual 
values  of  the  fields  from  these  equations.  These  require  the  knowledge  of  5q 
and  X,  which,  in  turn,  are  related  to  che  initial  values  of  the  rays  QP,  and 
P^P  as  well  as  che  attenuation  constant  *(j).  The  initial  value  of  the  field 
at  0 is  related  to  che  strength  of  che  source  by  L(0),  the  so-called  "launching 
coefficient,"  while  che  initial  value  of  the  field  at  is  related  to  the 
actual  field  on  che  surface  at  P^  through  the  "diffraction  coefficient"  D(P^). 

If  we  now  sum  up  che  contributions  of  all  che  modes,  we  obcain  the  final 
solution  [25]  for  che  field  radiated  in  che  shadow  region  by  an  infinitesimal 
magnetic  dipole  of  strength  >£  located  on  a smooth  convex  conducting  body 
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and  che  expression  for  G is  obtained  by  replacing  the  superscript  "h"  by 

"s"  in  (30b),  where  h and  s stand  for  hard  and  soft  boundary  conditions, 

h s h s 

viz.,  u«0  and  3u/3n  - 0,  respectively.  The  quantities  L ’ , D?’  and 

a*1’3,  appearing  in  (30b), in  general  depend  upon  the  local  geometry  and  the 
? 

electromagnetic  characteristic  of  the  surface,  frequency  k,  and  the  mode  of 
propagation.  They  are  determined  by  studying  the  asymptotic  expansions  of 
the  exact  solutions  of  some  special  canonical  problems.  Keller  and  Levy 
([20]  and  [27])  have  derived  the  first  few  terms  of  the  asymptotic  expansions 
for  D and  a by  considering  the  canonical  problems  of  scalar  diffraction  by  a 
circular  cylinder,  sphere,  elliptical  and  parabolic  cylinder.  A study  of  the 
above-mentioned  asymptotic  expansions  and  the  works  of  Franz  and  Klante  [28] 
and  Voltner  [29], who  have  also  investigated  the  same  problem,  as  well  as  a 
comparison  with  the  results  of  the  "direct  integral  equation  approach,"  reveals 
the  following  characteristics  of  the  solution:  i)  the  first-order  terms  in 
Che  asymptotic  expansion  of  D and  a are  independent  of  whether  the  problem 

i 

under  consideration  is  scalar  or  vector;  ii)  the  first-order  approximation  of 

f che  surface  along 

of  o , but  also  of 

a 

surface  transverse 
for  scalar  and 
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The  leading  terms  in  the  asymptotic  expansion  of  "diffraction  coefficient 
D,  "attentuation  constant"  a and  "launching  coefficient"  L are  presented  below: 
"Soft"  polarization: 
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'Hard"  polarization: 
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where  Ai(x)  is  the  Airy  function: 


Ai  (x) 


I f 


( cos(-|  + Xt)  dt  (x  real) 
o 


(31) 


(32) 


(33) 


(34) 


(35) 


(36) 


(37) 
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and  Ai(-r  ) - 0,  Ai'(-r')  ■ 0,  (Ai'  is  Che  derivative  of  Ai  with  respect 
? P 

to  its  argument).  Higher-order  terms  in  the  expansion  of  D,  a and  L have 
been  given  in  [25]  and  [30]  and  in  some  of  the  other  works  on  GTD  mentioned 
earlier. 

The  expression  (30)  is  convenient  to  use  in  the  lit  region.  In  the 
shadow  part  of  the  transition  region,  since  the  exponential  decay  of  the  terms 
in  (30)  is  weak.  The  convergence  of  the  series  representation  is  very  slow. 
Furthermore,  the  series  diverges  in  the  lit  part  of  the  transition  region. 
Consequently,  in  these  regions,  it  is  more  reasonable  to  use  an  integral 
representation  for  the  surface  ray  field,  which,  in  our  case,  can  be  expressed 
in  terms  of  Fock  functions  [25]. 

Attempts  have  been  made  to  establish  the  mathematical  validity  of  GTD 
and  to  minimize  its  "nondeductive  parts"  (parts  which  are  based  upon  physical 
intuition  or  the  study  of  the  asymptotic  solution  of  some  simple  problem 
geometrical  concepts  of  different  kinds  of  rays,  diffraction  coefficients, 
attenuation  constants,  etc.).  Kravtsov  [31]  and  Ludwig  [32]  have 
analyzed  the  field  near  the  caustic  surface  (smooth  envelope  of  a family 
of  ray) , and  have  developed  a "uniform  asymptotic  solution"  in  the  sense  that 
it  is  finite  at  the  caustic  and  reduces  to  geometrical  optics  away  from  the 
caustic. 


3.  Generalization  of  GTD  and  Investigation  of  Alternate  Methods 
3 . 1 Generalization  of  GTD  to  Arbitrary  Surfaces 

Keller's  generalization  of  GTD  for  the  analysis  of  the  field  diffracced 
from  a smooth  convex  object  is  closely  related  to  what  is  known  as  the 
"boundary  layer  technique"  in  the  theory  of  differencial  equations  [43].  On 
the  other  hand,  the  "uniform  asymptotic  theory"  is  analogous  to  the  method 
used  by  X.  E.  Danger  and  F.  J.  Oliver  to  find  the  asymptotic  solutions  of  the 


second-order  differencial  equations  near  their  "turning  points,”  which  are 
counterparts  of  the  transition  regions  in  our  case  [33],  [34],  and  [35]. 

The  procedure  is  based  upon  the  generalisation  of  the  geometric  optical 
interpretation  of  the  circular  cylinder  problem.  The  solution  obtained  by 
this  method  involves  some  functions  with  unknown  phase  and  amplitude,  similar 
to  Bessel  and  Hankel  functions.  Since  the  surface  of  a smooth  object  is 
actually  the  caustic  surface  of  diffracted  rays,  the  above-mentioned  formulation 
is  applicable  in  this  case,  too.  Lewis,  et  al.  [36]  have  modified  this 

solution  to  make  it  satisfy  the  boundary  condition  on  a convex  body.  Using 
ray  formalism,  they  have  obtained  an  asymptotic  solution  in  a complicated  form, 
which  they  call  "creeping  wave"  and  satisfies  the  boundary  condition  on  and  is 
uniformly  valid  near  and  away  from  the  surface.  It  should  be  mentioned  that 
the  method  has  been  developed  primarily  for  scalar  diffraction  problems. 

Creeping  waves  that  are  traveling  on  the  surface  of  the  body  generate 
other  kinds  of  diffracted  rays  in  the  presence  of  the  irregularities  in  the 
geometric  or  electromagnetic  characteristics  of  the  surface.  The  effects  of 
discontinuity  in  the  surface  curvature,  its  higher-order  derivatives, 
or  the  surface  impedance  have  been  studied  by  many  authors  [37],  [38], 

[39]  or  [40].  An  exhaustive  study  of  various  diffraction 

mechanisms  and  corresponding  diffraction  coefficients, 
and  constants  associated  with  the  propagation  of  creeping  waves,  has  been 
carried  out  by  Albertsen  [41]. 

At  this  point,  let  us  examine  one  of  the  most  important  features  of  GTD 
and  its  various  modifications.  GTD  formulation  is  essentially  scalar  in  nature 
and  is  heuristic  in  some  parts.  Thus,  when  GTD  is  applied  to  a vector  problem, 
it  is  not  surprising  that  the  coupling  between  various  components  of  the  fields 


are  neglected,  and  each  one  of  them  is  Created  as  an  uncoupled  scalar  wave. 

The  other  assumptions  in  GTD  are  concerned  with  the  directions  of  these  field 
components  and  the  kind  of  boundary  conditions  they  satisfy  (see  Sec.  2.4). 

As  mentioned  earlier,  non-deductive  parts  of  GTD  are  based  on  asymptotic 
expansions  of  known  solutions  to  some  selected  "canonical'  problems . Quite 
often  these  canonical  problems  are  not  general  enough  to  fully  and  accurately 
describe  the  local  behavior  of  the  field  for  an  arbitrary  structure.  Finally, 
most  of  the  canonical  problems  investigated  are  two-dimensional  in  nature.  The 
only  exception  to  this  is  the  sphere.  However,  in  30  far  as  the  geometric 
properties  of  the  surface  are  concerned,  the  sphere  is  a very  special  case 
since  its  radius  of  curvature  is  the  same  in  all  directions  and, consequently , 
the  surface  rays  are  torsionless.  Finally,  GTD  fails  when  the  observation 
point  is  located  in  the  transition  regions,  shadow  boundaries  or  in  the  neighbor- 
hood of  a caustic.  In  each  of  these  regions,  one  needs  to  carefully  modify  the 
GTD  formulas  and  often  such  a modification  is  not  too  simple.  Nevertheless, 
in  spite  of  these  difficulties,  GTD  is  recognized  to  be  a powerful  high-frequency 
technique  for  computing  the  leading  terms  of  the  asymptotic  solution.  Two  of 
the  principal  attributes  of  GTD  are  its  simplicity  and  wide  scope  of  application. 

3.2  Soectral  Domain  Approach 

We  now  examine  an  approach  different  from  GTD  which  uses  the  spectrum 

of  the  induced  current,  or  the  expression  for  the  radiated  field,  as  a starting 
point.  In  order  to  gain  a better  insight  into  the  curved-surface  radiation 

and  scattering  problem  and  to  verify  the  basic  assumption  of  GTD,  it  is 

worthwhile  to  consider  such  alternative  approaches,  particularly  if  they  apply 

to  canonical  problems  which  are  more  general  in  nature  than  those  employed  to 

derive  the  GTD  results.  An  example  of  such  a study  would  be  to  consider  the 


case  of  surface  ray  propagation  with  non-zero  torsion,  a situation  that  occurs 
when  a magnetic  dipole  source  radiates  from  a location  on  the  surface  of  a 
circular  cylinder. 

The  geomecry  of  the  problem  is  shown  in  Fig.  11  (p.  21).  The  radius  of  the 
cylinder  is  a and  the  source, which  is  an  infinitesimal  magnetic  dipole  with 
density  M,is  located  at  the  point  Q described  by  the  spherical  polar  coordinates 
(r«a,  $-0#).  Each  point  P on  the  surface  of  the  cylinder  is  defined  by  a 

"geodetical  polar  coordinate"  system  (a, 3),  where  a is  the  arclength  of  the 
geodesic  connecting  Q to  P and  3 is  the  angle  between  l (at  point  Q)  and 
geodesic  QP.  The  local  orthonormal  basis  vectors  Co, 3)  are  also  associated  with 
these  two  parameters.  The  observation  point  in  the  far  field  is  specified  by 
its  spherical  polar  coordinates  (r,3,$).  The  radiated  field  at  an  arbitrary 
point  can  be  expressed  in  terms  of  two  potentials , 5 and  % which,  in  cylindrical 
coordinates,  can  be  written  as: 
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For  the  problem  under  consideration,  we  can  express  the  spectral  weight 
coefficients  as 
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where 


In  order  to  derive  an  asymptotic  expansion  of  (37)  and  (38),  we  proceed 
as  follows.  As  a first  step,  we  apply  Watsonfe  transformation  to  the  infinite 
summation  with  respect  to  n and  employ  appropriate  asymptotic  formulas  for 
Hankel  functions  with  large  order  and  argument  to  derive  the  following  expressions 
for  (37)  and  (38)  under  the  conditions  that  ka  is  large  and  * small  compared  to  t; 


where 
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Fock's  functions  defined  in  Appendix  A. 


M and  M * components  of  M,  (M*n*0) 
4>  z 


!Iext,  applying  the  "saddle-noint"  technique  to  (42)  and  (43)  and  keening  only 

the  first-order  terms,  the  far  field  can  be  written  in  terms  of  its  components  along 

the  normal  and  tangent  to  the  surface  at  the  "stationary  point"  as 
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where 


P^:  is  che  scationary  point  of  which  turns  out  to  be  the 
same  as  che  point  of  diffraction  predicted  by  GTD. 


■lg  - (y)1/3  U-ir/2)  • sin173 


o ■ radius  of  curvature  of  geodesic  QP, 
o i 


J » arc  length  QP^ 


R » che  distance  between  the  point  of  diffraction  P^  and 
the  observation  point 


n^  * * 32;  normal  to  the  surface  at  P^ 


The  details  of  che  derivations  of  (44)  and  (45)  are  given  in  Appendix  3. 

Fig.  12  illustrates  the  geometric  meaning  of  some  of  che  parameters 
appearing  in  (44)  and  (45),  for  che  observation  point  is  located  in  the  shadow 
region.  In  this  case,  which  is  identical  to  l given  in  (16),  is  che 

reduced  distance  traveled  by  the  surface  ray  before  leaving  the  surface 
tangentially. 

In  the  lit  region,  che  geometric  interpretations  of  j and  % are  shown 
in  Fig.  13.  The  rays, like  OP^P,  that  do  not  obey  the  generalized  Fermat's 
principle  are  called  "psuedo-ravs"  [25].  The  ray  QP^P  appears  to  travel  along 
che  surface  up  to  che  point  P^  and  then  leaves  che  surface  at  P^  tangentially 
in  the  opposite  direction,  to  reach  the  observation  point  P.  It  should  be 
noted  that  formulas  (44)  and  (45)  give  us  the  contribution  of  che  ray  which 
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: Diffraction  of  rays  by  a cylindrical  oody. 


13:  Diffraction  of  "psuedo-rays . " 


travels  along  the  shortest  path  on  the  surface,  and  thus,  suffers  the  least 
attenuation.  It  is  not  difficult  to  see  that,  in  general,  there  are  infinitely 
many  rays  which  contribute  to  the  total  field  at  any  observation  point.  However, 
their  contributions  are  very  small  compared  to  that  of  QP^  and  their  phases  and 
amplitudes  can  be  determined  in  a similar  manner  (Appendix  3) . Several  other 
remarks  on  the  formulas  given  in  (44)  and  (45)  are  in  order: 

a)  Numerical  results  indicate  that  good  agreement  between  (49), 

(45),  and  the  exact  modal  solution  is  obtained  for  ka  > 10. 

b)  The  zeroth-order  terms  in  the  asymptotic  expansion  of  the  normal 
component  of  the  field  E are  identical  to  those  given  by  GTD; 

-1/3 

however,  the  k terms  derived  from  the  two  approaches  are 
different . 

c)  Tangential  component  of  the  field,  E,,  given  by  (45),  also 
is  different  from  the  corresponding  expression  based  on  GTD 
by  a multiplicative  polarization  factor.  Specifically, 


(M-$  ) 
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Consequently,  our  results  agree  in  GTD  only  for  the  circumferential 
ray,  i.e.,  for  3 » 0.  In  addition,  for  an  axial  r.agnecic  dipole 
(M  • a.  *0),  GTD  gives  a nonzero  value  for  the  field  in  the  3,  directic 
our  solution  predicts  that  this  field  is  identically  zero,  a result 
which  is  in  complete  agreement  with  the  exact  solution  for  the 


problem. 

In  contrast  to  GTD,  formulas  (44)  and  (45)  are  valid  in 
respect,  to  the  location  of  the  observation  point,  be  it  in  the 
lit,  shadow  or  transition  regions.  Although  not  valid  in  the  paraxial 
region  (3’90®),  they  can  be  generalized  to  work  along  this  direction. 
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Finally,  let  us  consider  the  possibility  of  the  generalization  of 
(44)  and  (45)  to  ocher  convex  surfaces  of  more  general  nature. 

3y  "more  general  surfaces,"  we  mean  chose  surfaces  which  are  not  substantially 
different  from  cylinders,  some  examples  being  cylinders  with  noncircular 
(ellipcical,  hyperbolic  and  parabolic)  cross  sections  and  conical 

surfaces  with  small  apex  angles.  The  key  step  in  a systematic  approach  to 

' 

generalizing  (44)  and  (45)  is  to  use  the  generalized  definition  of  5 given 
in  (16). 

Fig.  14  exhibits  some  initial  results  of  the  generalization  of  these 
formulas  to  the  case  of  a cone.  It  is  evident  that  results  obtained  from 
che  present  approach  agree  quite  well  with  che  series  solution  which  is  rather 
tedious  and  time-consuming.  We  also  observe  from  Fig.  14c  chat  there  is  a 
noticeable  discrepancy  between  the  analytical  solution  and  the  experiment.  Thus, 
within  the  range  of  experimental  error,  our  results  agree  quite  well  with 
those  published  in  the  literature. 

3 . 3 Approach  3ased  on  an  Asymptotic  Evaluation  of  the  Radiation  Integral 
of  the  Surface  Current 

As  a final  topic,  we  consider  an  approach  based  on  the  asymptotic 
evaluation  of  che  radiation  integral  expressed  in  terms  of  che  induced  surface 
current  which  is  itself  derived  in  an  asymptotic  manner  for  surfaces  with 
large  radius  of  curvature. 

It  was  shown  in  Sec.  2 chat  Fock's  theory  can  provide  us  with  an 
expression  for  the  scattered  field  in  che  neighborhood  of  a smooth  convex 
body  illuminated  by  a plane  wave.  Using  this  solution  in  conjunction  with 
the  reciprocity  principle,  we  can  find  the  far  field  radiated  by  a point  source 
located  on  che  surface  of  che  body.  3y  generalizing  the  definition  of  I in 
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. Measured  (Hughes) 
j Computed  (Hughes) 
,Computed(UI 


Theta 


Comparison  between  spectral  domain  results  (UI) , 
modal  approach  (Hughes) , and  experimental  measure- 
ments (Hughes) . The  UI  results  are  derived  from  a 
generalized  version  of  (44)  and  (45)  for  a cone. 
The  Hughes'  results  have  been  reproduced  from  [63] 
and  are  based  on  a modal  series  of  13  terms.  All 
results  are  for  \/2  radial  slot  on  a cone  of 
half-angle  10s. 
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Fock's  theory,  we  can  also  write  the  final  result  in  a GTD  format  and  represent 
it  as  a surface  ray.  The  total  field  at  a point  on  the  surface  is  obtained 
by  adding  all  the  possible  rays  which  reach  the  observation  point  ?.  Various 
techniques  can  be  used  to  determine  the  field  propagation  along  these  rays. 

For  instance,  when  the  source  is  located  on  the  surface,  and  the  surface  is  a 
conical  one,  the  field  at  each  point  can  be  decomposed  into  two  parts. 

F - T + F2  (47) 

where  F^  is  the  geometrical  optics  field  when  the  observation  point  is  directly 
illuminated  by  the  source,  and  is  the  creeping-wave  contribution  derivable 
via  an  extension  of  Fock's  theory  when  the  point  is  in  the  shadow  region.  The 
other  term,  F^,  is  the  so-called  tip  contribution,  and  can  be  obtained  by 
physical  optics  or  GTD.  Goodrich  et  al.  [42]  have  applied  this  procedure  to 
find  the  radiation  pattern  of  slot  arrays  on  cones. 

The  approximate  induced  surface  current  distribution  cam  be  obtained  by 
Fock's  theory,  GTD  [13],  [14],  [16]  and  [25]  or  some  other  appropriate  high 
frequency  technique.  The  induced  surface  currenc  due  to  a magnetic  dipole  on 
a perfectly  conducting  circular  cylinder  and  cone  has  been  calculated  by  Chang 
et  al.  [44],  and  Chan  et.  al  [45]  whose  procedure  is  based  upon  an  asymptotic 
expansion  of  the  exact  modal  solution  to  the  above-mentioned  problems.  Lee, 
et  al.  [46]  and  [22]  have  treated  the  same  problem  by  a method  based  on 
Fock's  asymptotic  solution  of  the  problem  of  a sphere  [47].  These  expressions 
for  the  current  distribution  can  be  used  in  the  radiation  integral  representation 
of  the  far  field. 

The  numerical  evaluation  of  this  integral  is  a formidable  task,  especially 
when  the  frequency  is  very  high.  Thus,  it  is  highly  desirable  to  have  an 
analytical  and  explicit  formula  for  the  far  field  expressed  in  terms  of  the 
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surface  current.  We  now  discuss  an  approach  for  accomplishing  this  task  and 
examine  the  problem  of  deriving  an  asymptotic  expansion  of  the  far  field 
radiated  due  to  a point  source  located  on  the  surface  of  a smooth,  conducting, 
and  convex  body  of  an  arbitrary  shape. 

Consider  an  arbitrary  smooth  convex  surface  S shown  in  Fig.  15.  Let 
a magnetic  dipole  source  be  locaced  at  a point  Q on  S.  We  parametrize  the 
surface  S introducing  a "geodetical  polar  coordinate”  system  with  the  pole 
located  at  Q such  that  an  arbitrary  point  P^  on  the  surface  is  defined  by  a 
pair  of  numbers  (a, 3),  where  is  the  arc  length  of  the  geodesic  QP^  and 


is  the  angle  between  QP^  and  some  reference  direction  at  Q.  Unit  vectors 


along  the  constant  parameter  curves  a and  3 are  locally  orthogonal.  The 
unit  normal  to  the  surface,  n,  is  given  by  n - o * 3.  An  element  of  length 
in  this  coordinate  system  may  be  written  as 


2 2 2 
ds  “do  + G dS 


(48) 


The  radiation  integral  for  the  scattered  far  field  can  be  written 


iuu 


E * “vjr1  J 1 
s 


. «*ElzHgl  ds 


(49) 


where  R is  the  distance  between  any  point  on  the  surface  and  the  observation 
point.  In  the  geodetical  polar  coordinate  system,  we  can  rewrite  a scalar 
component  of  (49),  3ay  M,  in  terms  of  a double  integral  of  the  following 
general  form 


:<  - / ; rw.a.P)  .0  4,  a 


(50) 


J 
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where  we  have  assumed  Che  following  form  for  Che  surface  currenc: 


J(PL)  - tCP^  exp(-jkLj)  - J^j  + J33 

■ k - jk?;  k - 2ir/\,  k2  > 0,  k2  <<  k 


(51) 


where  f(P^)  is  a relacively  slowly  varying  funccion  when  k,  che  free -space 

wave  number,  is  large.  This  assumpcion  is  based  upon  a close  scruciny  of 

differenc  asympcocic  formulas  given  for  che  induced  surface  currenc. 

Then  che  observacion  poinc  P is  locaced  in  che  shadow  reeion, 

main  concribucion  co  (51)  comes  from  a small  neighborhood  of  che  scacionary 

poinc  of  che  incegrand ,and  che  scacionary  phase  mechod  for  mulciple  incegrals 

([48]  and  [49])  is  applicable.  The  asympcocic  expansion  of  (49)  has  been 

-5/3 

derived  up  co  che  order  k . The  decails  of  calculadon  have  been  presenced 
in  Appendix  C.  The  final  resulc  is 


d - e"jW4  • e076  • r ci/2)  • r ci/3)  • o 2 
u 0 


6776  • r(2/3)  • r (1/2)  . 0 473 

J 


A,  B,  and  C are  dependent  upon  geometric  properties  of  the  surface  at  the 
stationary  point  which  turns  out  to  be  exactly  the  same  as  the  "point  of 
diffraction"  of  surface  rays.  The  quantities  A,  B,  and  C are  given  by 


, . 0 o r3B 

4 m -i-Q  _ —2.  . — (— ) + —2.  r_ 

A 2 G 3a  2 3o  ' 2G  lc> 

0 0 


2 

+ (L30)2  - (l/2>^4]  + 0 
3o~ 


where 


B - L3<7/G172,  C - - l/o. 


o ■ radius  of  curvature  of  the  geodesic 
o 

°g  * geodetic  radius  of  curvature 

„33  „ 3o  A e. 


L3  , L”J  ■ coefficients  of  the  second  fundamental  form  of  the  surface  (S) 


A geometric  interpretation  of  these  parameters  has  been  illustrated 

1/2 

in  Fig.  16.  It  i3  evident  from  this  figure  that  ) 1 is  simply  the 

divergence  factor  of  the  rays  leaving  the  surface  tangentially  at  the  point 
of  diffraction.  In  using  formula  (36) , we  should  bear  in  mind  that  the 
various  terms  in  and  (J^  are  not  of  the  3ame  order.  For  example,  in  the 
deep  shadow,  is  exponentially  larger  than  Ja. 


The  formulas  given  in  (56)  have  been  tested  and  compared  with  ocher 
available  solutions.  An  important  conclusion  derived  from  this  comparison 


is  that  although  the  method  of  radiation  integral  is  based  on  less  restrictive 
assumptions,  it  is  perhaps  not  as  useful  as  the  spectral  domain  approach 
because  the  stationary  point  of  the  phase  of  the  integrand  in  (50)  is  of  the 
second  order,  and  hence,  the  asymptotic  expansion  of  this  integral  converges 
rather  slowly  except  when  kp^.  is  very  large  ( = 40  or  more). 
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APPENDIX  A:  FOCK  FUNCTIONS 


In  studies  of  radio-wave  propagation  around  the  earth  by  Van  der  Pol, 
3remmer,  Pryce,  Fock,  and  others,  and  also  the  later  studies  of  diffraction 
of  electromagnetic  waves  by  certain  bodies  of  revolution  ( [51] , [ 52 1 , [53 ] , 

[54] , [551 , [48] , [56] , [57] , [58] , [59] , [60]  and  [50]),  a class  of  universal 
functions  was  introduced  which  can  be  used  to  predict  the  amplitude  and  the 
phase  of  the  reflected  or  diffracted  field  by  smooth  convex  surfaces  [17]. 

An  exhaustive  treatment  of  these  functions  which,  in  general,  are  defined 
as  Fourier  integrals  having  combinations  of  Airy  integrals  in  their  inte- 
grands, has  been  carried  out  by  N.  A.  Logan  [61].  (See  also  3ovman,  et  al. 

[1]  and  Logan  and  Yee  [17]). 

Since  the  first  extensive  application  of  these  functions  to  diffraction 
theory  was  done  by  Fock,  many  auchors  named  them  after  him.  Here  we  list 
only  the  most  important  formulas  and  expressions  for  these  functions  without 
going  through  the  details  of  their  derivations.  We  have  followed  Logan's 
set  of  notations  for  these  functions  [61].  However,  since  his  time  dependence 
factor,  exp(-iuit),  is  different  from  one  we  have  used  throughout  this  paper, 
namely  exp (■‘•jwt) , our  expressions,  listed  below,  are  conjugates  of  what 
have  been  presented  in  [ 61] . 

First  we  3tart  with  general  definitions.  Fock's  most  general  form  of 
the  "Van  der  ?ol-Bremmer  diffraction  formula"  is 


V(x,yi>y2,q)  - exp ( j*/ 4) 


jv(t  - v<) 


v' (t)  - qv(t) 
w 2 ( t ) - qw,(t) 


W,(t  - y>)  ; dt 

i 


(A.  1) 


where  w^(t) , w?(t),  u(c)  and  v(t)  are  rock. -type  Airy  functions,  defined 


u(t)  * v t 3i(t) 


, v(t)  - /rr  Ai(c) 


w^Ct)  - u(c)  + jv(t)  , w2(t)  * w^(t) 

We  note  that  w^  and  can  also  be  defined  as  in  Sec.  2.  y>  and 

y<  are  the  larger  and  smaller  of  the  two  numbers  y^  and  y^.  V(x,y1ty7,q)  is 
proportional  to  the  attenuation  suffered  by  an  electromagnetic  wave  generated 
by  a source  located  at  reduced  height  y^  above  the  surface  of  a smooth  con- 
vex body,  whan  it  reaches  the  observation  point  located  at  reduced  height 
y^  above  the  same  surface,  x is  the  reduced  distance  between  the  source  and 
the  observation  point  along  the  surface,  and  q is  dependent  upon  the  impedance 
of  the  surface.  Let  us  consider  some  useful  limiting  cases. 


When  y^  * y7  * 0,  then  V(x,0,0,q)  is  denoted  by  V^,  where 


Vg(x.q) 


jir/4  , — r » e " xC  w,(t) 

/ x 1 2 
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7(t)  - qw7(t) 


(A.  2) 


We  also  have 


VQ(x,0) 


( 2 

lim  i -2jxq  VQ(x,q) 
- 


rjxc  w2(t) 


■*  3/2  r 

— • x 


(A.  3) 


When  y^  » 0 and  y,  ■*  »,  then  7 -*•  V1  (x,q)  : 


» -j  xt 
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V1U,,>  * 3?  ' i-  ■-UjCO 


and  also 


?(x)  - V^(x,0)  * — 
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(A.  7) 


3ased  on  Equations  (A. 6)  and  (A. 7),  a class  of  functions  can  be 
defined: 


f(n)00 


;M(x) 


* T 


f ta  • e^  . a dnf(x) 

ir  *2(c) 


(A.  3) 


, .i'\n  f „n  *jxt  jn  / \ 
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wi(t) 
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dx 


where  7 is  any  path  in  the  complex  t-plane  which  comes  from  -®  in  a sector 
defined  by  - t _<  arg(t)  < - y and  goes  to  +°°  in  the  sector  -y  < arg(t)  < y. 

In  what  follows,  we  will  give  the  suitable  formulas  for  f(x)  and  g(x)  in 
different  ranges.  Tabulated  values  and  graphs  of  these  functions  can  be 
found  in  [ 57  J , f 54 ] and  [61]. 

When  x is  very  large  and  negative,  the  following  asymptotic  expansions 
for  f(x)  and  g(x)  can  be  used  [61]: 


4 x3/3  f j 1 4 

f(x)  - -2:xejX  1 + -±r  + ~ - J- 


175 


395  j_318l75 
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4x  lx''  64x'  16x  ~ I024x 

g(x)  - 9e^3/3/l  - -3 L + 1£°2  + J™L . 11122 

■ \ , 3 6 ,.9  ,12 

1 +X  X 04X  04X 

The  above  formulas  are  valid  and  accurate  for  x <<  -1.  For  moderate 
values  of  x,  namely,  -1  <_  x <_  1,  it  is  difficult  to  find  an  appropriate 
expression.  Although  there  are  some  analytical  techniques  like  "stationary 
phase  method"  or  "Poisson  summation  formula"  which  may  be  used  to  evaluate 
f ' and  g for  these  values,  another  possible  way  which  is  probably  easier 
and  more  efficient  is  to  interpolate  the  tabulated  values  of  these  functions 
in  this  range. 
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In  che  vicinicy  of  zero  (jx|  =0),  che  Taylor  expansion  can  be  used 
co  calculate  f and  g.  The  coefficients  are  given  by 


f (n)  (0)  - a~j  (5mr/6-ir/3)  . ^ 
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where  t(\,u)  is  che  generalized  "tau"  function: 
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:(X.u)  - l — -v  , X > 1 
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(49n2  + 364n  + 39349/16) / (29  • 32) 


When  x is  large,  and  positive,  residue  series  can  be  used  to  compute 
:(n)  and  g<n)  : 


f ^ (x)  * e 
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APPENDIX  B:  DERIVATION  OF  FORMULAS  (44)  and  (45) 


Here,  we  consider  only  Che  derivation  of  the  asymptotic  expansion  of 
■}  for  a circumferential  magnetic  dipole.  In  this  case,  i may  be  written 


where 


r a 

(2ir) 2 J- 


S(kt)  « l e 


-jk  z 

a 2 • S(kc) 
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-to  . H”)('v> 
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Applying  the  Watson  transformation  to  (B.2), 


(■  II  (k  a)  -jv(p-ir) 

Q Tic  ^ * ■**“  • "• — 1 • ■ ■ ■ . ■ ■ • 

C 2 J'  C+D  H(2)(k  a)  Sin  V1T 
v t 


(3.3) 


where  C and  D are  shown  in  Figure  17.  Or, 
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Substituting  the  expansion 
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where  p^  * p and  p£  * 2*  - p,  in  (B.4),  the  result  will  be: 


2 ; (•  » H;  '(k  p)  -jv(P,+2*l) 
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Each  term  of  the  above  expansion  is  associated  with  a "creeping  wave" 
travelling  in  a counterclockwise  (i  * 1)  or  clockwise  (i  * 2)  direction 
around  the  cylinder.  Following  the  ray  concept,  each  creeping  wave 
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Incroducing  a new  incegracion  variable  a: 


we  have: 
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R^  *<z2  + [p  + a^  + 2*1  - ir/2)]|'1^2 


Now  (B.10)  cakes  the  following  form: 
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(B. 15) 


Y is  Che  pach  of  incegracion  in  Che  complex  a-plane,  which  is  shown  in 
Figure  18. 

Now  we  deform  Che  pach  of  incegracion  inco  che  "sceepesc  descenc  pach," 
SDP,  passing  Chrough  Che  saddle  poine  of  Che  phase  of  Che  inCegrand.  Per- 
forming che  "saddle” poinc  incegracion,"  we  can  derive  che  asympcocic 
expansion  of  (B.15)  for  large  kR^.  The  firsc  order  cerm  is: 
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where  R_, , and  - . , are  che  values  of  chese  oaramecers  ac  che  scacionary  poinc 
ils  ils  • r 

specified  by  a - 3^, • 
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Re  (a) 


Eqn.  (B.16) is  Che  creeping-wave  representation  of  the  far  field.  If  the 
cylinder  is  large  (lea  » 1)  and  j *|  is  not  very  close  to  it,  then  only  the 
first  term  (i  - 0,  1-1)  has  the  most  important  contribution  to  the  total 
infinite  sum,  and  the  ocher  terms  are  not  significant.  Neglecting  the  other 
terms,  we  obtain  the  result  given  in  (44  and  45).  It  should  be  emphasized 
that  (44)  and  (45)  are  not  valid  when  ]3j  is  close  to  t/2  (paraxial  region), 
because  in  this  case,  k£a  is  very  small,  and  (B.7),  (B.S)  and  (3.9)  no  longer 
apply. 

The  other  formulas  can  be  derived  in  a similar  manner. 
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APPENDIX  C:  ASYMPTOTIC  EVALUATION  OF  THE  RADIATION  I1TEGRAL 


Consider  the  following  double  integral: 


UOO  - f [ g(x,y)  • dxdy 

JDJ 


(C.l) 


where  g(x,y)  is  rather  slowly  varying,  and  *(x,y)  has  a stationary  point 


(x^.y^)  inside  domain  D.  The  objective  is  to  derive  an  asymptotic  expansion 


for  (C.l)  when  k is  large. 


Suppose  g and  0 have  the  following  forms  around  (x  ,y  ): 
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N.  Chako  [<*8]  has  derived  the  following  asymptotic  series  for  U: 
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. 


m+n>l 


ar*ti>l 


la  order  Co  apply  chls  procedure  co  che  incegrals  of  the  type  (5.4)  for 


which 


"•Then  F is  one  of  the  components  of  J(1  - RR) , one  should  first  determine 


the  stationary  point  of  ft,  wherein  its  first-order  derivatives  vanish.  The 


;hi3  point,  and  then  insert  them  into  (C.3).  We  just  give  che  main  formulas 


polar  coordinate  system.  As  discussed  previously, in  this  system,  z is  the 


chen  ve  have  the  following  sec  of  relations 


3323a 


_G_  3G/3o  32G/3c2 

R2  R 2 


333  33G/33  3G  3G  _1_  . 3a  33  1 32G 

3g3  " 2R  33  ‘ 3c  ’ 4G  L L 2 3a38 


(C. 16) 


(C. 17) 


where  is  the  radius  of  curvature  of  the  geodesic. 

Equations  (C.12)  and  (C.13)  determine  the  location  of  the  stationary 

point.  Ac  this  point  R ■ x^,  which,  if  we  introduce  the  ray  concept. 

cells  us  chat  the  surface  rays  leave  the  surface  at  Che  "point  of  diffraction" 

tangentially.  Equation  (C.14)  indicates  that  the  stationary  point  is  of 

second  order,  so  chat  we  need  higher-order  derivtcies  of  the  phase.  L , L" 

33 

and  L are  coefficients  of  the  second  fundamental  form  of  the  surface 
evaluated  at  the  stationary  point.  They  are  defined  as 

L ' So  ' *3  • L ' \e  ' *3  ' L ’ *36  ' *3 

Using  the  relationships  given  above,  one  can  find  the  expansion  co- 
efficients g,  . , a , b and  A in  (C.3).  Zeroth  and  first-order  terms  in 
Scl  mn  mn  pq 

(C.3)  give  us  formulas  (5.6). 

A few  remarks  should  be  made  concerning  the  expansion  presented  in 
(C.3).  First  of  all,  (C.3)  is  a doubly  infinite  series;  therefore,  for  each 
fixed  power  of  k 2 a finite  number  of  terms  should  be  summed  up.  The 
coefficients  of  various  terms  in  these  finite  sums,  namelv  A 's  , become 

pq 

very  complicated  when  p and  q are  greater  chan  0 or  1.  Another  difficulty  with 
this  series  is  chat  when  the  stationary  point  of  the  phase  is  of  an  order 
higher  than  1,  the  difference  between  the  order  of  the  successive  terms 
(when  they  are  ordered  according  to  the  descending  power  of  k)  becomes  very 
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small,  and  consequently  Che  infinite  series  converges  very  slowly. 


For 


instance,  in  cur  problem  where  5=3  and  t ■ 2 (stationary  point  is  of 
second  order),  sometimes  the  difference  becween  the  orders  of  successive 
terms  is  , which  indicates  the  weak  convergence  (in  an  asymptotic 

sense)  of  the  expansion  in  the  cases  where  the  frequency  is  not  very  large. 
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